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DIFFERENTIAL ABSORPTION UDAR MEASUREMENTS OF ATMOSPHERIC WATER VAPOR 
USING A PSEUDONOISE CODE MODULATED AIGaAs LASER 

BY 

Jonathan A. R. Rail 
ABSTRACT 

Lidar measurements using pseudonoise code modulated AIGaAs lasers are 
reported. Horizontal path lidar measurements were made at night to terrestrial targets 
at ranges of 5 and 13 km with 35 mW of average power and integration times of one 
second. Cloud and aerosol lidar measurements were made to thin cirrus clouds at 13 km 
altitude with Rayleigh (molecular) backscatter evident up to 9 km. Average transmitter 
power was 35 mW and measurement integration time was 20 minutes. An AIGaAs laser 
was used to characterize spectral properties of water vapor absorption lines at 
811.617, 816.024, and 815.769 nm in a multipass absorption cell using derivative 
spectroscopy techniques. Frequency locking of an AIGaAs laser to a water vapor 
absorption line was achieved with a laser center frequency stability measured to better 
than one-fifth of the water vapor Doppler linewidth over several minutes. Differential 
absorption lidar measurements of atmospheric water vapor were made in both integrated 
path and range-resolved modes using an externally modulated AIGaAs laser. Mean water 
vapor number density was estimated from both integrated path and range-resolved DIAL 
measurements and agreed with measured humidity values to within 6.5% and 20%, 
respectively. Error sources were identified and their effects on estimates of water 
vapor number density calculated. 
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1. INTRODUCTION & BACKGROUND 


Water vapor comprises <3.0% of the Earth’s atmosphere but is extremely 
important to both life processes and atmospheric physics. The importance to atmo- 
spheric physics is defined in a strategic research plan 1 developed during an October 
1990 workshop held by the Global Energy and Water Cycle Experiment (GEWEX). 

Water vapor is: (1) by virtue of its latent heat transfer property, a principal medium 
by which energy is exchanged among the components of the Earth system, i.e. atmo- 
sphere, hydrosphere, cryosphere, and biosphere; (2) the predominant greenhouse gas, 
playing a crucial role in many radiative processes which regulate the global climate; and 
(3) essential in many atmospheric processes, e.g. cloud formation and precipitation, 
which determine climate variations, especially on regional scales. Currently, our lack 
of knowledge of the distribution of atmospheric water vapor and its variability prevents 
reliable assessment of potential regional or global climate change. This deficiency may 
be corrected through a combined campaign of observations and modeling of atmospheric 
water vapor. 


1.1 Water vapor & global circulation modeling 

Global circulation models (GCM’s) attempt to parameterize all significant 
processes which drive atmospheric circulation. Since water vapor plays a major role in 
most of these processes, parameterization of moist processes is critical. Understanding 
these complex water vapor processes and developing parameters which accurately 
describe them requires a detailed knowledge of atmospheric water vapor distribution and 
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its variability. The desired spatial resolution is horizontal grids less than 200 km on a 
side and 1 km vertical layers up to the tropopause. The measurements should be 
temporally resolved to permit sensing of diurnal variations. Total column-content 
measurements are also needed to initialize these models with realistic total atmospheric 
water vapor content. 

1 .2 Current water vapor measurement techniques 

A global upper air balloon sounding network provides regular radiosonde 
measurements of atmospheric water vapor.2 These radiosondes have numerous short- 
comings. Geographical spacing of sounding sites is typically 400 km in North America 
and western Europe and sparse elsewhere. There is little or no radiosonde coverage over 
the oceans. Temporal resolution is also poor with measurements made at twelve hour 
intervals. Variation in sensor types and reporting practices used in different countries 
leads to significant uncertainties in humidity data. Almost no data is collected in the 
upper troposphere since measurement accuracy degrades at temperatures below -30°C 
and most sondes do not report in conditions less than -40°C.^ 

Satellite observations of water vapor provide more complete geographical 
coverage, improved temporal resolution, better horizontal resolution, and in some cases 
better vertical resolution compared to radiosondes. Current operational satellite 
techniques include passive infrared and microwave radiometers in both vertical and 
limb-scanning modes. 

Passive infrared radiometers typically use several broadband infrared channels 
in the 6-14 p.m region to observe water vapor radiance's from different layers of the 
atmosphere. An example is the High Resolution Infrared Spectrometer (HIRS/2) 
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instrument of the Tiros Operational Vertical Sounder (TOVS) investigation flown on the 
TIROS-N and NOAA 6 through 10 satellites. 4 The water vapor content at three vertical 
levels of the atmosphere may be calculated from these radiances. Unfortunately radiance 
observations are inhibited by cloud cover, therefore, water vapor profiles cannot be 
retrieved from overcast areas and are only estimated over partially cloudy areas. When 
used in a limb-scanning mode, radiometers such as the Stratospheric Aerosol & Gas 
Experiment (SAGE II) instrument on the Earth Radiation Budget Satellite (ERBS) can 
provide 1 km vertical resolution profiles of water vapor in the upper troposphere and 
stratosphere. 5 However, due to the tangential path geometry of the measurement, the 
horizontal resolution in this mode is poor. Also, since the limb-scanning or solar- 
occultation mode requires specific alignment of the sun, earth, and spacecraft the 
measurements have a limited coverage. 

Passive microwave sounding instruments, e.g. Scanning Multichannel Microwave 
Radiometer (SMMR) on the Nimbus-7 satellite, make total column measurements of 
atmospheric water vapor over the oceans by measuring the differential absorption 
strength of water vapor between 18 and 21 GHz.6 The accuracy of the column content 
measurements is considered to be better than ± 1 0%. 1 However, vertical profiles of 
water vapor cannot be obtained and measurements are not made over land, sea ice, or 
areas of precipitation. 


1 .3 Laser remote sensing of water vapor 

There has been significant sustained interest in using laser remote sensing 
techniques to make measurements of atmospheric water vapor at high spatial and 
temporal resolution. The ultimate goal is making spaceborne measurements of 
atmospheric water vapor profiles and total column content. Most laser remote sensing 
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measurements of atmospheric water vapor have used either ground based Raman or 
differential absorption lidar (DIAL) approaches. 

1.3.1 Raman lidar 

Measurements of atmospheric water vapor profiles using Raman lidar were first 
reported by Melfi, Lawrence, and McCormick 7 and Cooney 8 using a second harmonic of a 
ruby laser and by Strauch, Derr, and Cupp 9 using a N 2 laser. The Raman scattering 
technique involves detecting radiation which has been shifted in wavelength due to 
interaction with the scattering molecule. This wavelength shift, or Stokes shift, is equal 
in energy to a vibrational-rotational or rotational transition in the scattering molecule. 
The backscattered power of the wavelength shifted signal is proportional to the 
concentration of Raman scattering molecules and inversely proportional to Thus, 
the primary advantage of Raman lidar is that it offers a direct measure of species 
concentration or mixing ratio by comparing the Raman signal of water vapor to the 
Raman signal of N 2 or 0 2 . However, Raman scattering is a very weak process and the 
Raman backscatter signal for water vapor is typically two to four orders of magnitude 
weaker than the elastic backscattered signal. This creates problems with isolating the 
inelastic scattering signal from the elastic backscatter. 10 Also, the weak scattering 
cross section typically limits water vapor Raman lidar to nighttime measurements at 
ranges of less than 10 km. To increase the Raman signal and make daytime measure- 
ments, high power lidar systems have been developed to operate at wavelengths from 
248.5 to 268.5 nm. Unfortunately at these wavelengths absorption by molecular oxygen 
and ozone can attenuate the transmitted beam and solar irradiance can obscure the 
backscattered water vapor Raman signal. 
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1 .3.2 Water vapor absorption spectroscopy 


The water vapor DIAL technique employs elastic scattering of radiation by 
aerosols and molecules and absorption by water molecules to estimate the number 
density of atmospheric water vapor. The near-infrared spectrum of water vapor has 
many absorption lines with absorption line strengths spanning several orders of 
magnitude. There are three prominent water vapor absorption bands in the near 
infrared, centered at 720, 820, and 940 nm. Figure 1 is a plot of atmospheric 
transmission of solar irradiance calculated using the LOWTRAN database 11 and the 1976 
U.S. Standard Atmosphere. 


WAVENUMBER (cmA-1 ) 

13500 1Z600 11700 10800 



Figure 1 . Atmospheric transmittance of solar irradiance from 0.7 - 1.0 jim. Vertical 
path through entire atmosphere as calculated with LOWTRAN 
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The strong absorption feature at 760 nm is due to the oxygen (0 2 ) A-band. The water 
vapor bands centered at 720, 820, and 940 nm correspond to electronic ground state, 
vibration-rotation (VR) transitions of the water molecule, a semi-rigid rotator, 
asymmetric top structure belonging to the symmetry point group c 2v . 12 The 720 and 
940 nm bands have been extensively studied 13 ’ 14 - 1 5(1 6,17,18,1 9 w hile the 820 nm 
band has received less attention 20 The prominent VR transitions in the 820 nm region 
originate in the v^yj = 211 vibration state and terminate in the vibration ground state, 
V 1 Y 2 Y) = 000. Figure 2, calculated using the HITRAN database and the 1976 U.S. Standard 
Atmosphere, shows atmospheric transmission for a wavenumber range 12000-12400 
c/rr 1 over a 5 km horizontal path at sea-level. 

WAVELENGTH (um) 


0.830 0.825 0.820 0.815 0.810 

— i 1 1 1 1 



Figure 2. Atmospheric transmission spectra vs wavenumber showing individual water 
vapor lines for a 5 km horizontal path at sea level. 
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Figure 2 shows that individual water vapor absorption lines in this region are resolved. 
The specific line assignments are determined by the upper and lower state rotational 
quantum numbers J' ,K' a , K' b & J" ,K" a , K" b , respectively, and are shown in Figure 3. 



Figure 3. Line assignment and relative intensities of strong water vapor lines in the 
(2,1,1) band. 21 

The line assignments have been abbreviated using the notation (/’ T . —J" T «) where 
T’=(K' a -K' b ). The R, £),and P branches correspond to A/=1,0,-1, respectively. 

The absorption cross section a for a particular water vapor absorption line is a 
function of wavelength, pressure, and temperature and can be expressed as the ratio of 
the line strength S(T ) to the Lorentz linewidth y(p,T) 22 
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oa, P ,T) = Si P~~ • /a - A m >, 

n-y(pj) 


(1 ) 


where /(A - A m ) is the wavelength dependence of the cross section referenced to 
maximum absorption at A m . If the laser linewidth is one-fifth or less the width of the 
absorption line, then /(A - A m ) may be taken as unity. In DIAL measurements, the on- 
line laser should be spectrally narrower than the absorption line of the atmospheric 
species of interest to maximize measurement sensitivity. Ponsardin et al. reported that 
the laser linewidth must be no larger than one-quarter of the absorption linewidth to 
ensure that no systematic error is made when estimating water vapor number density 
from DIAL lidar data. 23 Water vapor absorption lines in the 820 nm region typically 
have pressure broadened linewidths, at STP, of 1-5 GHz full width at half maximum 
(FWHM). 

Atmospheric temperature and pressure vary with altitude, therefore the 
dependence of the absorption cross section on temperature and pressure as well as 
pressure shift of the absorption line center can be a significant source of error in water 
vapor DIAL measurements. Pressure broadening effects can be dealt with by using 
spectrally narrow lasers while pressure shift effects can be minimized by locking the 
laser to line center at lower pressures. This compensates for pressure induced shifts of 
the line center by locking the laser at a vapor pressure equivalent to a higher altitude. 
Temperature sensitive parameters have been derived for the 720 & 940 nm water 
vapor absorption bands. 1 4 ’ 1 6>1 81 9 

The 820 nm water vapor absorption band is attractive for atmospheric DIAL 
measurements due to the abundance of both strong and weak water vapor lines and the 
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availability of both photon counting detectors and high power AIGaAs laser diodes in this 
wavelength region. A large dynamic range of line strengths is required for DIAL 
measurements to accommodate various measurement ranges and varying water vapor 
concentrations. 


1.3.3 Differential absorption lidar 

Schotland first proposed measuring vertical distributions of atmospheric water 
vapor using DIAL techniques in 1966. 24 The DIAL technique employs two wavelengths to 
estimate atmospheric water vapor number density. One wavelength X on is selected to 
coincide with the center of a water vapor absorption line while a second wavelength X a ff 
is selected to fall in a nearby nonabsorbing region. If X on and X Q jj- are within a few 
cm 1 of one another, then the elastic scattering properties of the atmosphere are 
assumed to be identical and can be neglected. Laser power at both wavelengths is 
transmitted into the atmosphere (either simultaneously or sequentially) and is 
elastically scattered by molecules and aerosols into the field of view (FOV) of the lidar 
receiver. Provided that the transmitted laser power is low and the receiver FOV is 
small, multiple scattering effects may be neglected. The received power P(X l ,R) of the 
scattered light, at each wavelength, X on and X off , and range R, is given by the single 
scattering lidar equation 25 


P&L.R) 


- P ^O 


P L 


R 


2 ^sys 


Z(R)P n (X L ,R)AR e~ 2 l° K( ^’ R)dR . 


( 2 ) 


Here P L is the transmitted laser power, A a is the area of the telescope, R is the range, 
x sys is the lidar system optical transmission, fi n (X L ,R) is the atmospheric backscatter 
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coefficient at range R and X L , £(/?) is the overlap function of the laser divergence with 
the receiver FOV, A R = cr L /2 is the length of the atmospheric measurement cell, x L is 
the receiver integration time, and k(X l ,R) is the total attenuation coefficient due to 
aerosol and molecular scattering and absorption. 



Figure 4. Overlap of laser transmitter and receiver FOV, (A) zero overlap, (B) 
partial overlap, and (C) complete overlap. 

Biaxial lidar systems have separate transmitter and receiver optical paths which 
require that the transmitter and receiver be carefully aligned to one another. The 
overlap function, £(/?), gives the range dependence of transmitter energy crossing into 
the receiver FOV. The overlap function of a biaxial lidar is graphically depicted in 
Figure 4. 

Equation (2) may be rearranged and simplified to give the expected received 
lidar signal in photoelectrons per second for a homogeneous atmospheric path. For most 
measurement ranges of interest, there is complete overlap (^ = 1). The total 
attenuation coefficient for a homogeneous path can be broken into its constituents and the 
expected on-line photoelectron count rate N on becomes 


( 3 ) 


^pA « 


Here 77 is the quantum efficiency of the detector, E ph = hc/X L is the laser photon 
energy, k m is the water vapor absorption coefficient, <J m is the molecular scattering 
coefficient, k a is the aerosol absorption coefficient, and o a is the aerosol scattering 
coefficient. Mean values for the volume backscatter coefficient, fi n (X L ,R), and 
individual attenuation coefficients have been measured and are catalogued as functions of 
wavelength, altitude, latitude, and season and are available from several sources. 26,27 

The DIAL technique can be used to estimate the number density of water vapor 
molecules at a specific range (range resolved) or the average for a path. The water 
vapor absorption cross section, cr(cm 2 ), and water vapor absorption coefficient, 
k m (cm -1 ), are related by 25 


k m {X L ,R) = n{R)-o{X L ,R\ (4 

where n(R) is the number density of water vapor molecules in the atmosphere at range 
R. 


When the laser is tuned off of the water vapor absorption line, the water vapor 
absorption coefficient k m = 0, and the expected photoelectron count rate is 


N 0 ff{h off ,R)- 




E ph R 


2 %ys Pji 'R)^R ■ e 


-2 R{o m +k a +o a ) 


1 1 


( 5 ) 


The backscatter coefficient, P n (X L >R) t which is assumed to be constant between X ^ 
and X 0 ff , can be estimated from lidar measurements and the visibility or meteorological 
range representing the atmospheric conditions at the time of the measurement. The 
average number density of water vapor molecules, n(AR ), in a single range bin defined 
as AR = R 2 - R v can be expressed in a simplified DIAL equation. The DIAL equation is 
derived by taking In (N on /N off ) for a range R , and subtracting ln(N on /N off ) for an 
adjacent range R 2 , and solving for the water vapor absorption coefficient k m . Using the 
difference in absorption cross section Acr, the mean number density of water vapor 
between /?, and R 2 can be expressed as 


n(AR) = 


1 

2A(J'M 


Wjjg 2) 


( 6 ) 


Here A<7= 0(X on )- a{X off ) is the difference in water vapor absorption cross section 
between X on and X off , N on = N on T and N off = N off T are the integrated photoelectron 
counts at X on and X off respectively, T is the lidar measurements integration time, and 
/?, and R 2 are ranges which define the boundaries of A R, the range bin of interest. 

Simulations of spaceborne water vapor DIAL measurements at night indicate that 
accuracy’s of <10% with horizontal resolutions of 100 km and vertical resolutions of 
less than 1 km can be achieved with a 5 Hz, 150 mJ pulsed laser at 727 nm.^ A nadir- 
pointing DIAL instrument using an alexandrite laser transmitter has recently measured 
water vapor profiles from a NASA research aircraft. ^9 



Currently operational water vapor DIAL and Raman lidar systems use solid state, 
gas, dye, or excimer laser systems. 30 These laser transmitters are large, bulky and 
inefficient with complex optical pumping systems. Due to these limitations, these laser 
transmitters are not viable candidates for remote sensing applications where size, mass, 
and power are constrained. 


1 .4 AIGaAs lasers 

Semiconductor AIGaAs laser diodes are much smaller and more efficient laser 
sources than there solid-state counter parts. Current single element, single 
longitudinal mode lasers are capable of reliably producing 200 mW of continuous power. 
Recent developments of master oscillator-power amplifier (MOPA) devices have boosted 
this CW power to the 1-10 W level. With continued improvements in power levels, 
MOPA devices based on AIGaAs semiconductor lasers, can be considered for aircraft and 
spaceborne DIAL remote sensing instruments. The emission wavelength of AIGaAs 
devices, between 780 and 860 nm, is well matched both to atmospheric water vapor 
absorption bands and to photomultipliers and silicon detectors. 

Some of the earliest operational lasers were gallium arsenide (GaAs) homo- 
junction diodes which were cryogenically cooled. 3 ^ However, due to high threshold 
current densities, cooling requirements, and erratic lifetimes, these early laser diodes 
were not practical for many applications. The introduction of aluminum to these 
devices, creating AI( X )Ga(-|. X )As compounds and the development of the double 
heterostructure (DH) geometry resolved both the high current density and reliability 
problems which plagued these early lasers. 


AIGaAs has the useful property that the bandgap energy can be varied over a wide 
range by changing the AI( X )Ga(i- X ) ratio, with only a negligible change in lattice para- 
meter. The lattice parameter, a 0 , changes only 0.14 % when changing AI( X )Ga ( i. X )As 
from AlAs to GaAs. This substantially improved the reliability of AIGaAs lasers by 
reducing lattice mismatch while permitting a broad manufacturing range of the emission 
wavelength. 

A potential barrier or quantum well which simultaneously confined the injected 
carriers and created a rectangular optical waveguide was created by burying the double 
heterostructure laser in undoped AIGaAs material and tailoring the bandgap. The 
waveguide confined the laser emission and reduced absorption losses in the GaAs 
substrate. Subsequent improvements in epitaxial layer growth, wafer material purity, 
and lattice matching resulted in lower threshold current densities and subsequently 
longer lived lasers. Current commercially available quantum well AIGaAs lasers are the 
smallest, most efficient, and least expensive lasers available, and have demonstrated 
lifetimes of greater than 30,000 hours. 

However, single mode laser diodes are limited in the peak optical power that they 
can produce due to the possibility of catastrophic facet damage. Catastrophic facet 
damage is mechanical failure or melting of the facet due to intense optical fields. ^ ^ The 
width of the emitting region and the pulse length are important factors which determine 
the power density level (watts/cm 2 ) at which failure occurs. With this limited 
capability to produce high peak powers, conventional short pulse lidar measurement 
techniques are not practical using single mode AIGaAs lasers. These lasers are, however, 
well suited for pseudonoise (PN) code modulation, which has a nearly 50% duty cycle 


and whose noise-like correlation properties permit range resolved backscatter 
measurements. 32 


1 .5 Pseudonoise (PN) codes 

Maximal length pseudonoise codes (PN codes) are the longest non-repeating 
series of ones and zeroes that can be generated by a digital shift register of a given 
length. 32 Digital shift registers with feedback can be used to produce maximum length 
PN sequences of 2"- 1 bits, where n is the number of stages of the shift register, and a 
bit is a single element (1 or 0) of the sequence. Long shift register code generators, 
typically 8 to 12 stages, produce more useful code lengths, from 255 to 4095 bits. PN 
code modulation affords range resolved measurements with range resolution determined 
by the modulation rate and a maximum unambiguous range determined by the code length. 

A simplified block diagram of a PN code lidar system is shown in Figure 5 and is 
described as follows, a PN code modulated laser beam is transmitted through the 
atmosphere and a small fraction, less than one photon per bit, is backscattered into the 
receiver. The receiver detects these photons and synchronously accumulates a 
photoelectron count over 10^ to 10® repetitions of the code. This received histogram is 
a record of counts versus time delay (range bins). The histogram represents the 
convolution of the atmospheric backscatter function with the transmitted code. 

Correlating the histogram with a stored version of the PN code yields the atmospheric 
lidar signal. 

In PN code lidar, the correlation function is used to compute the lidar signal from 
the detected backscattered photons. The noise-like correlation properties of the PN code 
permit recovery of the lidar signal by cross-correlating the received histogram with 


the transmitted sequence. The received stream of photons has been modified by the 
atmospheric path and target. Hence, the correlation function, which computes the lidar 
backscatter signal vs. range, contains atmospheric backscatter and absorption 
information as a function of range. 


Atmospheric aerosol & molecular 
scattering 



Figure 5. Simplified block diagram of PN code lidar system. 

1 .6 PN code lidar background 

Several previous lidar systems have used PN code modulation to obtain range 
resolved signals with low peak power lasers. In 1983, Takeuchi et al. externally 
modulated an Argon laser, at 514.5 nm, with a PN code to measure the lidar return from 
a smoke plume at 1 km. 33 in 1986, Takeuchi et al. demonstrated a PN code aerosol lidar 
using a single 30 mW AIGaAs laser diode transmitter. 34 This system measured lidar 
returns from falling snow, smoke, nighttime aerosols, and cloud structure. In 1988, 
Norman and Gardner proposed a PN code technique for performing laser ranging 
measurements to satellites, and presented a signal and error analysis. 3 ^ In 1992, 













Abshire et al. reported nighttime measurements to tree canopies using a PN code 
modulated AIGaAs laser. 36 With a modified version of this system, Rail et al. reported 
measuring nighttime aerosol profiles to 4.0 km altitude and cloud returns to 8.6 km 
altitude. 37 In 1993, Abshire and Rail reported a simplified PN code lidar theory and 
nighttime measurements to cirrus clouds at 13.5 km and terrestrial targets at 13 km. 38 

Chapter 2 follows with a thorough analysis of PN code lidar theory and operation. 
In addition, the performance of an AIGaAs water vapor DIAL system is predicted based on 
the single scattering lidar equation and the system parameters of an existing breadboard 
lidar. 


Chapter 3 reviews water vapor absorption spectroscopy. Criteria for selecting 
near-ideal absorption lines for water vapor DIAL measurements are developed. 
Candidate absorption lines are selected and profiled with a laser diode in a multipass 
absorption cell using derivative spectroscopy. Spectroscopic parameters including 
linewidth and line center wavelength are measured for the candidate absorption lines. In 
addition, a technique for actively stabilizing the laser diode frequency to a water vapor 
absorption line is developed. 

Chapter 4 describes two prototype AIGaAs lidar systems that have been assembled 
and tested. The first system described is a cloud and aerosol lidar and the second is a 
wavelength stabilized, water vapor DIAL system. System diagrams and instrument 
characterization are included for each lidar. 

Chapter 5 presents lidar data acquired with each system and subsequent data 
analysis. Lidar measurements over horizontal paths to terrestrial targets are compared 
with theory using the single scattering lidar equation. Water vapor density is estimated 


from DIAL measurements and is compared with ground based humidity measurements 
made at local airports. Frequency stability of the actively stabilized laser diode is 
estimated. 

Chapter 6 summarizes results of the laser diode frequency stabilization effort, 
the cloud and aerosol lidar experiments, and the water vapor DIAL measurements. 
Future work on AIGaAs lidar and altimetry is discussed. 


2. PHOTON COUNTING AIGaAsUDAR THEORY 


The theory governing the operation of AIGaAs lidar may be developed by 
considering a laser transmitter which is intensity modulated with a maximal length 
pseudonoise (PN) code, of length m, and the single scattering lidar equation. An 
expression for the detected lidar signal may be derived by assuming a single reflector or 
scatterer at a fixed range. This expression may then be generalized by extending it to 
multiple scatterers at distributed ranges. The lidar, Figure 5, transmits the PN code 
modulated light into the atmosphere and a fraction of these photons are scattered into the 
FOV of the lidar receiver and detected. The detected photons produce a sequence of 
photoelectron emissions which, for a single scatterer, occur in the range bins of the 
original PN code but lagging the original code by a time delay corresponding to the 
roundtrip range delay. These photoelectron emissions are accumulated into a histogram 
and stored in memory. Cross-correlating the histogram with the original PN code yields 
the atmospheric lidar signal. 

2.1 PN codes and their properties 

PN codes have three noise-like properties which are important to their use in 
ranging and lidar measurements^ 9 : (1) The number of ones and zeroes are nearly equal, 
always one more one than zero. (2) The distribution of ones and zeroes in a sequence is 
well defined and always the same from one sequence to the next. (3) The normalized 
autocorrelation function of the sequence yields unity correlation for zero relative delay 
and near zero correlation for all other values of delay. 

1 9 



A nearly equal number of ones and zeroes allows the transmitting laser to operate 
with a -50% duty cycle, with the peak operating power twice the average power. Laser 
diodes are well suited to operating at such high duty cycles. The distribution of ones and 
zeroes within a PN code sequence determines its noise-like correlation and spectral 
properties. Although maximal length PN sequences do repeat (and are deterministic), a 
sampling of ones and zeroes within the sequence is nearly random and can be made 
arbitrarily close to random by simply increasing the sequence length. The auto- 
correlation function of a PN code measures the degree of agreement between a code and a 
time delayed replica of itself. 

A maximal length PN code, a,, has elements ^=(0,1), where i = 0,...,m-l. The 
sequence may be generated by an n-stage shift register where the code length m is 
related to n by 

m = 2 n - 1 . (7) 

An example of a 7-bit PN code, generated by a three stage shift register, is shown in 
Figure 6. 



Figure 6. Three stage shift register and the 7-bit code it produces when the initial 
state is the all ones state. The feedback taps are added modulo-2 with an EX-OR gate. 
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An m-bit code sequence has (m + 1)/2 ones and ( m - 1)/2 zeroes. The m- b i t 
code, a it with elements (0,1) may be expressed alternatively as the code, a\, with 
elements (-1, 1) generated by 


a i = 2a i ~ 1 • (8 

This alternate form of the code is useful in lidar where the amplitude of the correlation 
function contains information regarding atmospheric backscatter and extinction 
properties. An example of a 255 bit, amplitude shifted PN code is shown in Figure 7. 



Figure 7. Example of a 255 bit a\ PN code. 



2.1.1 Autocorrelation & cross-correlation 


The autocorrelation function measures the degree of agreement between a code 
sequence and a time delayed replica of itself. The autocorrelation function, of a 

function, x[«] , is defined as 40 


+oo 

$ xx (^) = * 

m -- °© 


(9) 


Computing the autocorrelation function for the code, <x* v yields^ 


m- 1 




i '=0 


f 1 J = o 

1/m y*0. 


( 10 ) 


where y is modulo- tn. In PN code lidar, the cross-correlation function is used to 
compute the lidar signal from the detected backscattered photons. The cross-correlation 
function, <pxy[n], of the functions x[n] and y[n] is defined as 40 


+ « 

$xy( n ) = Xm+n • 
m=-~ 


(ID 


It is important to note that the cross-correlation function possesses odd symmetry such 
that 


$yx [ • 


( 12 ) 
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Computing the cross-correlation function for the code sequences rj,- and o.\ yields the 
correlation function ^.(y) which is shown in Figure 8 


, / , \ , \(m + l)/2 j = 0 


(13) 





Figure 8. The correlation function ^-(y) of a PN sequence of length m. 


The peak value of the correlation function occurs at zero time shift to the original code, 
j = 0, and its amplitude is equal to the number of 1‘s in the code, (rw + l)/2. For other 
values of time delay y , the correlation function is zero. The cross-correlation function, 
Qaa'U)’ is more easily visualized in Figure 9 with a 7-bit PN code. 
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zero time delay (j=0 ) 



<= iOi H2\3 \4 \5 \ 6\ 7 \ 8\ 

'L a i a ' i+ j = 1+1+1 +0+0+1 +0 = 4 = (m+1)/2 

1 bit time delay (j=-l) 



'La.a\ ml = 1 + 1 + 0 + 0 - 1 +0 -1 =0 


wrap around last 
bit 



/= I Q I / I 2 I .? I 4 I 5 I 6 I 7 I 8 I 


wrap around last 
n bits 


'La i a’ i _ n = 1 + 0 + 0 - 1 + 0 + 1 - 1 ... = 0 


Figure 9. Sample cross-correlation function, <f> aa ’{j), for a 7-bit PN code. Note that 
a\ has been amplitude shifted as in Eq. 8 and that for time delays other than zero, the 
“extra" n bits of a\ correlate with the first n bits of a r 
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2.1 .2 Maximum range and range resolution 


If the transmitter modulation rate or bit frequency is denoted f b , then the bit 
period or time duration for each bit is t b = 1 j f b . The lidar system range resolution is 

A* = (c/ 2),,, (1 

where the factor of two is due to the roundtrip traveled by the light. The code sequence 
length m and the bit period t b determine the maximum unambiguous range which may 
be measured with an m-bit code sequence. The maximum range is given by 


x = m-(c/2)-t b . (15 

Since the sequence repeats identically after /n-bits, if is too small, there can be 
an ambiguous situation where light is scattered from two ranges, separated by one-half 
of the code sequence length, ( m -AR)/2 . This constitutes a “wrap-around” of 
consecutive code sequences and the return (correlation peak) from the more distant 
target occurs in the same range bin as the closer target. To avoid this, m must be 
selected which has a maximum unambiguous range, greater than the anticipated 
maximum range. 

2.2 Expected AIGaAs lidar signal & signal-to-noise ratio 

If the laser is operating with an average power, P Q , then the transmitted laser 
power of the i th bit of a, is 
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The expected received power P , scattered from range Rj, is governed by the single 
scattering lidar equation and can be written as 


P(A,L,Rj) = P 0 --$; • x sys ■ £(R)Pn(A-L’Rj)&R ' e ° ’ ( 1 7 ) 

R j 

where P a is the transmitter average laser power, A 0 is the effective area of the 
receiver telescope, Rj = j ■ A R is the distance to the scatterer, is the lidar receiver 
optical transmission, £(R) is the laser divergence and receiver FOV overlap function, 
P x (X l ,Rj) is the atmospheric backscatter coefficient {km~ x • $r _1 ) at wavelength X L and 
range Rj, A/? = (c/2)- r b is the range cell size, and is the total attenuation 

coefficient ( krrr 1 ) due to scattering and absorption by aerosols and molecules. The 
expected received power may be converted into an expected photoelectron count rate 
using the energy of the transmitted photons and the quantum efficiency of the detector. 
The instantaneous photoelectron count rate, N jt produced by a photon counting receiver 
is 

-2jK(X L ,R)dR 

Zph “j 

where N denotes the instantaneous photoelectron count rate for range bin j, p,_j is the 
instantaneous transmitted power of the /'* bit reflected from the j th range bin, i] is the 


quantum efficiency of the detector, = hc/X L is the laser photon energy, and the 
overlap function £(/?) has been assumed to be unity. The expectation value of the 
photoelectron count rate at the bit of the PN code scattered from the j A range cell 
therefore may be written as 



= p H ■ y 


R ? 


-2 jK(A L ,R)dR 

■ARe 0 


( 19 ) 


where 


7 = 


V-A 0 r : 


sys 


he, 


lc A 


( 20 ) 


is the combined lidar system parameters. Equation 19 may be simplified by defining 


h(X L ,Rj) -2/rtA t .*W 
= — - — s — —ARe 0 


R 


( 21 ) 


as an atmospheric scattering and extinction function. The expression for the expected 
photoelectron count rate from the j * range bin, N it becomes 

NiJ=P H rGj. (2 

In addition to the signal photoelectrons, a term for background light must be included. 
This term, b it is the background photoelectron count rate in the i* bin due to 
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unmodulated laser transmitter light and solar radiation scattered into the FOV of the 
receiver. This background photoelectron count rate can be expressed as 


bi = J V • t b dX, (23) 

where S b (X) represents the spectral radiance of the sky background, AX is the spectral 
width of the receiver’s band pass filter, is the receiver transmission efficiency, 

£2 0 is the receiver acceptance solid angle, A 0 is the area of the receiver telescope 
aperture, and t b is the bit period. The total expected photoelectron count rate in the i‘ h 
bin due to signal photoelectrons scattered from the j‘ h range cell and background 
photoelectrons is 


Nij = Pi-j-r-Gj + b t . (24) 

These photoelectron counts are accumulated into receiver bins synchronously with the 
transmitted code sequence, creating a histogram of received counts over the integration 
period. Each receiver bin corresponds to one bit of the PN code sequence. If the signal is 
accumulated over L cycles of the code sequence, then the receiver integrates for L • t b 
seconds at each bin. The total integration time is T = m L t b seconds. The total 
integrated counts in the i A bin is then 

Nij = Lt b p H YGj +L-t b b i . (25) 
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This gives the functional form of the lidar return for a single scatterer. The extension to 
multiple scatterers at distributed ranges is done by summing over all range bins j, 
where j is modulo - m yielding 


m - 1 


N i =Lt b Y^Pi_jGj+Lt b - bi 
7=0 


(26) 


N i represents the histogram of received counts in each of i receiver bins due to signal 
scattering from range cells y = 0. ..m — 1 and background counts. The histogram contains 
information about the atmospheric path or transfer function which transformed the 
input photons to the received photoelectrons. The lidar signal, S n , may be generated by 
cross-correlating the histogram, N it with the modulation code, a\. 


m - 1 

S n = 'Z N i a 'i~n Y 

1 = 0 


m-lm-l m - 1 

E lLPi-j G j a 'i-n + Y,k a 'i- 
i=0 j = o i=0 


(27) 


where n is modulo- m and the property of odd symmetry Eq (12) has been used. 
Substituting pi_j = 2 P 0 ■ a^j yields 

m-l/n-l m- 1 

S n = 2 P 0 L • t b ■ y J X Oi-jGjCf^ +L t b ■ • 

i=0j=0 i=0 

Exchanging the order of summation, 
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( 29 ) 


S n = 2 P 0 L t b ■ yX 'L a i-j G j a ’i~n +L t f 'LW-n . 
y=oj=o i=o 


and using the properties of PN code cross-correlation 41 from Equation 13 


<t>aa'U- n ) = 


m - 1 


i = 0 


[(m + l)/2 

I 0 




(30) 


yields. 


m-1 _ 

S n =2PoL<b-r'Z<P aa 'U-n) G j +b - (31) 

j = o 

The background count rate, b x , has been assumed independent of the modulation code a- 
and therefore replaced with a time averaged background count, b = L-t b -bi which can be 
dealt with independently of the signal photoelectrons. The cross-correlation is non zero 
for n = j only, yielding the lidar signal 

S n = P 0 Lt b y(m + l)G n . (32) 

Equation 32 represents the lidar signal (without integrated background counts) as a 
function of range bin index n which is identical to range index j used earlier. The 
response function or atmospheric response may be derived by solving this equation for 
G n once S n has been calculated from the histogrammed data. 
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Due to the nature of light detection, the photoelectron counts in the received 
histogram Eq. (26) have a Poisson distribution. The probability function of the 
Poisson distribution is 


Ppi-X’H) — 



( 33 ) 


The average count x is an estimate of the mean, \i , of the Poisson distribution. The 
mean is found by calculating the expectation value of x 




x^—e * 
x =o\ * ! y 


,X-\ 


■"‘■'Si 


(34) 


where the last summation in Eq (34) is the series expansion of the exponential 



(35) 


The counts in each bin of the histogram Eq (26) have a Poisson distribution. The 
photoelectron counts fluctuate from observation to observation simply because random 
samples of events distributed randomly in time contain numbers of events which 
fluctuate from sample to sample. 42 For Poisson distributions, the mean equals the 
variance and the standard deviation is the square root of the variance. To compute the 
lidar signal for each range bin, Eq (29) sums the counts from each histogram bin over 
the entire received sequence. Therefore the variance of each lidar signal range bin 
(correlation space) is the sum of the variances of each histogram bin (histogram space) 
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which is the sum of all signal and background counts accumulated. Consequently, all 
range bins of the lidar signal S n Eq (32) have the same variance given by 


m - 1 

var(Sn) = o?. = 2 JV i- (3 6) 

i-0 

The signal-to-noise ratio (SNR) of the lidar signal in any range bin n is the signal 
count in that range bin (minus the average background count) divided by the standard 
deviation, a, 


SNR(n) = 



- y(m + l)G„ 

m- 1 

i = o 


(37) 


2.3 Performance calculations 

The expected performance of a single-color AIGaAs lidar is estimated using the 
single scattering lidar equation Eq. (18). These performance estimates are used to 
predict the performance of a differential absorption water vapor AIGaAs lidar. The lidar 
system parameters used were based on the existing prototype AIGaAs lidar and are listed 
in Table 1 . The atmospheric conditions, including backscatter coefficient and total 
extinction coefficient, were adapted from the 1962 U.S. Standard Atmosphere and the 
AFGL MidLatitude-Summer model. The total extinction coefficient, k(X l ,R), can be 
expressed as the sum 

/c( , R ) — k m + o m + + (T a , ( 38 ) 
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where k m is the water vapor absorption coefficient, a m is the molecular scattering 
coefficient, k a is the aerosol absorption coefficient, and or, is the aerosol scattering 
coefficient. Values for the individual attenuation coefficients have been adapted from the 
AFGL MidLatitude-Summer model based on the U.S. Standard Atmosphere (USSA) 1962 
and are catalogued as functions of wavelength, altitude, latitude, and season. 26 The 
individual attenuation coefficient values for 0.86 pm are listed in Table 2 as a function 
of altitude. 


Table 1 

AIGaAs lidar prototype system parameters. 


Laser power 35 mW Average 

Laser wavelength 810 nm 

Photon energy @ 810 nm 2.44E-1 9 J/photon 

Transmitter divergence 100 prad 

Telescope Area (20 cm diam./8 cm obscur.) 0.026 m 2 

Telescope Field of view (FOV) 160 prad 

Detector detection probability 0.25 PE/photon @ 830 nm 


Table 2 

Calculated off-line values for individual attenuation coefficients at 0.86 pm as a function 
of altitude, latitude, and season. 


Molecular (H20) 

Aerosol 

(clear) 

Height (km) 

k m 

(km~') 

o m [km- x ) 

k a (km ~' ) 

o a (*m _1 ) 

0 

< 

10E-6 

1.93E-3 

1.52E-2 

9.03E-2 

1 

< 

10E-6 

1.84E-3 

1.01E-2 

5.99E-2 

2 

< 

10E-6 

1.66E-3 

4.41E-3 

2.61E-2 

5 

< 

10E-6 

1.22E-3 

5.58E-4 

3.30E-3 

1 0 

< 

10E-6 

7.11E-4 

3.10E-4 

1.83E-3 


Note: Values for attenuation coefficients are for Midlatitude Summer conditions 


33 















The original USSA 1962 did not contain a water vapor distribution which was added 
later. 43 The water vapor number density used to calculate the water vapor attenuation 
coefficient 26 is listed in Table 3. 


Table 3 

Atmospheric water vapor density vertical profile. 


Height (km) 

Pressure 

(mbar) 

Temperature 

(K) 

Water vapor 
(q/m 3 ) 

Number density 
(molecules/cm 3 ) 

0 

1013 

294 

14.0 

4.60e1 7 

1 

902 

290 

9.3 

3.1 Oel 7 

2 

802 

285 

5.9 

1 ,96e1 7 

5 

554 

267 

1 .0 

3.33e1 6 

1 0 

281 

235 

0.064 

2.1 3e1 5 


The on-line water vapor attenuation coefficient k m is calculated for an absorption line at 
811.617 nm, adjusted for pressure and is listed in Table 4. 


Table 4 

Atmospheric water vapor attenuation coefficient k m vs altitude. 


Height (km) 

i 

E 

0 

0.44 

1 

0.33 

2 

0.24 

5 

0.059 

1 0 

0.007 


The received photoelectron count rate was calculated for single color horizontal 
path lidar measurements at altitudes of sea-level, 1, 2, 5, and 10 km and one-way path 
lengths of 1, 2, 5, and 10 km. Estimated values for the volume backscattering 
coefficient 44 and the off-line total attenuation coefficient are listed in Table 5. When 
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the laser is wavelength tuned to the center of an absorption line, the molecular 
absorption term, k m , dominates the total attenuation coefficient, k(X l ,R). However, 
the other attenuation coefficients may not be neglected. As an example, the absorption 
line at 811.617 nm was used to estimate photoelectron count rates for horizontal path 
DIAL measurements made at the center of the absorption line. The line strength for this 
line is 2.54 xlO -24 cm~\ / molecule ■ cm~ 2 which makes it a moderately strong absorption 
line in the 820 nm band. The total on-line attenuation coefficient was calculated using 
the water vapor absorption coefficients listed in Table 4 and are listed in the last column 
of Table 5. 


Table 5 

Backscatter coefficient and total attenuation coefficients j < 0 g & vs altitude. 


alt(km) 


K off (km l ) 

K on (km ! ) 

0 

2E-6 

0.1074 

0.547 

1 

4E-7 

0.0718 

0.401 

2 

2E-7 

0.0321 

0.272 

5 

IE-7 

0.0050 

0.064 

1 0 

6E-8 

0.0028 

0.010 


The estimated off-line photoelectron count rates vs altitude and horizontal range are 
listed in Table 6. These estimated photoelectron count rates are applicable as the off- 
line DIAL measurement of tropospheric water vapor. The estimated photoelectron count 
rates vs altitude and horizontal range for the on-line channel are listed in Table 7. 
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Table 6 


Expected off-line photoelectron count rates for various length horizontal path aerosol 
lidar measurements vs altitudes. 


Off-line 

Altitude (km) 

Horizontal Range (km) 

1 

2 

5 

1 0 

0 

1.12E + 05 

2.26E + 04 

1.90E+03 

1.62E+02 

1 

2.41 E + 04 

5.22E+03 

5.43E+02 

6.62E + 01 

2 

1 .31 E + 04 

3.06E+03 

4.04E+02 

7.31E + 01 

5 

6.89E+03 

1.71E + 03 

2.65E+02 

6.29E+01 

1 0 

4.15E + 03 

1 .03E + 03 

1.62E+02 

3.94E + 01 


Table 7 


Expected on-line photoelectron count rates for 811.617 nm water vapor absorption line 
for various length horizontal path lidar measurements vs altitude. 


On-line 

Altitude (km) 

Horizontal Range (km) 

1 

2 

5 

1 0 

0 

4.52E + 04 

3.66E+03 

2.00E+01 

1 .80E-02 

1 

1 .22E + 04 

1 .34E + 03 

1.82E+01 

7.40E-02 

2 

8.03E + 03 

1 .1 6E + 03 

3.57E+01 

5.71 E-01 

5 

6.12E + 03 

1 .34E + 03 

1 .46 E+02 

1 .91 E + 01 

1 0 

4.09E + 03 

1 .00E + 03 

1.51 E+02 

3.39E+01 


An estimate of the relative measurement error inherent to the AIGaAs DIAL 
technique can be derived by looking at the statistical fluctuations of the photoelectron 
count about the mean. Taking a ratio of the expressions for N on and N off , Equations (3) 
and (5), we find 


N„ 


N 


exp 


off 


-2jk m {R)dR\, 


(39) 
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where the contents within the brackets is the optical thickness due to attenuation by 
water vapor absorption. Rearranging this equation yields a relation between water 
vapor absorption and the on and off line photoelectron counts 


R 

\k m (R)dR = 


0 


'T ln 


N. 


on 


N. 


off 


( 40 ) 


The statistics of the water vapor absorption integral depend upon the statistics of the 
ratio of N 0n /N 0 ff. We can describe these fluctuations for N on and N a ^ as 


Non=(N m )' 




( N °») 


(41 ) 


and 


N «r -{"*)■ 


1 + 



( 42 ) 


The equations may be simplified by defining the fluctuating part as 


_ Won 


and 


e - ^off 

ff " K) ’ 


( 43 ) 


( 44 ) 


respectively. Since e on and e Q ff fluctuate symmetrically about the mean, 
( £ on) = { £ off) = 0- For large values of N, 
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( 45 ) 


AN 

(N) 


«!<=>£«!. 


Therefore the variance of the fluctuating part, e, is 

var < £ )=77 


and the standard deviation is 



Forming the ratio 

Non _ (Q-( 1 + £ J 

N off ( N off )( 1 + £ off) 

and simplifying, 


(46) 


(47) 


(48) 


Non _ (Q-Q + ggn) y 

N off ( N off ) ■ (l + £ off ) (l “ £ 0 ff) 


yields an expression for N 0n /N 0 ^ 


N. 


on _ 


<*«.>•(> 


N 


off 


+ £ on Zqff £ on-off J 

{ N off)'( l+e2 °ff) 


(50) 


Since e on and e off «1, the cross product and square terms may be ignored leaving 
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The mean value (n o „/N oj ^ may be formed yielding 



( 52 ) 


because ( e on ) - (toff) - 0- The variance of ( N 0n /N 0 ff ) is the square of the expectation 
value {^on/^off) multiplied by the sum of the variances of N on and N c 


var 


No 

\ N °ff . 


' (NS 

vK), 


(var(e < ,„) + varfe^)) 


(53) 


Substituting expressions for variances from Equation 46 yields, 


var 


'jO ({ N on 


N. 


Off ) 


vKr), 


' 1 X ^ 

1 

'y N on N Q ff J 


(54) 


The standard deviation is the square root of the variance, 


a 

v 


N n 


N. 


off J 


■■ (W«.) / 1 | 1 

K) 


(55) 


and the relative error in the ratio N 0n /N 0 ff is defined as the standard deviation divided 
by the mean, 


( 56 ) 


i^on) | 

1 1 4 1 - 

<T.M 

{ N on ) ( N off) 
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which yields 
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For ease of plotting, assign / = N 0n /N 0 jj , where / represents the water vapor 
absorption, then the expression for relative error becomes 
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rel. error = 



(58) 


This expression for relative error of N 0n /N 0 ff vs. N OJ gr. the off-line photoelectron 
count rate, is plotted in Figure 10 for 0.1 < / <1.0. Figure 10 indicates that for DIAL 
measurements with less than 10% relative measurement error and off-line photo- 
electron count rate greater than 1 000 counts/sec the fractional absorption / must not 
be less than 0.1. That is, greater absorption of the on-line channel by water vapor will 
lead to a larger relative error. However, the figure is misleading. It suggests that 
weaker absorption lines, which will yield less fractional absorption, will also yield a 
smaller relative error. 
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Figure 10. Plot of relative error of N 0n /N 0 ff vs. N for 0.1</<1.0. 


But as the fractional absorption approaches 1 , 
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A more important error to be determined is the relative error in the estimate of 
the water vapor number density, n(AR). The functional form for this error may be 
derived by using Eq (6), the DIAL equation, and the expression for the standard deviation 
of the ratio N 0n jN 0 ff, Eq. (55). The DIAL equation is 


n(AR) = 
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The desired expression is the relative error of n(AR). Eq (60) has the form: 

/j(A/?) = Cln[x] . (61) 

where C is a constant including the prefactors in Eq (60), and x = uv. Here 
u = N on {Ri)/^off { R \)< and v = N off( R i)/ N on i R 2)- The quantities u and v are assumed 
independent of each other. The relative error of n(AR) may be expressed as 



where 
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and 

— = C — ln(jc) = C - . 
dx dx x 


(64) 


Substituting the expressions x = uv, Eq (63) & (64) into Eq (62) yields the 
expression for the standard deviation of n(AR) 


a n = C 
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The relative error is formed by taking the ratio of Eq (65) with Eq (60), the 
expression for n(AR) yielding 
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( 66 ) 
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This expression may be simplified using the properties of the natural logarithm function 
yielding 
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For small range bin sizes, N on ( 1) = N on ( 2) and N ~ N 0 ff(2) and Eq (67) becomes 


-£l = 

n 2 In / ’ < DO > 

where / = N 0n jN 0 fj , the fractional absorption due to water vapor. The expression in Eq 
(68) is plotted in Figure 11. As the absorption due to water vapor decreases, i.e. / => 1 
the relative error blows up and Eq (68) reduces to 
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Conversely, as the absorption due to water vapor increases, f =$ 0, Eq (68) becomes 
indeterminate. Both the numerator and denominator of Eq (68) approach infinity. 
Applying L'Hopital’s rule to Eq (68) and simplifying yields 



n 



(70) 


which blows up to infinity as /=> 0. 



Figure 11 . Relative error of water vapor number density estimate as a function of off- 
line photoelectron counts and fractional absorption due to water vapor, / = N 0n /N 0 jf . 

Also plotted in Figure 11 are four data points from the estimated performance of 
the prototype AIGaAs lidar system. The on-line and off-line photoelectron count rates, 
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Tables 6 & 7, are based on system parameters listed in Table 1 and the 811.617 nm 
absorption line. The “X" indicates the expected lidar performance at 5 km altitude and 5 
km horizontal range; the marks the lidar performance at 5 km altitude and 2 km 
horizontal range; the “A” indicates lidar performance at 2 km altitude and 2 km 
horizontal range; and the “0” marks the expected lidar performance at sea-level and 2 
km horizontal range. The 811.617 nm line is too strong to make low altitude (<2 km) 
DIAL measurements of water vapor over horizontal ranges greater than 2 km. However, 
the line is too weak to make short path DIAL measurements above 2 km altitude, where 
water vapor density is lower and / => 1. Also, at higher altitudes where the backscatter 
coefficient is much weaker, the average laser power is too low to make long range (>5 
km) DIAL measurements. For a water vapor DIAL lidar to be accurate (minimal relative 
error) and have a practical range of up to several kilometers, Figure 1 1 indicates that a 
more powerful laser transmitter (higher N off ) and both weak and strong absorption 
lines are required. However, too strong of an absorption line will completely attenuate 
the on-line signal and the relative error will blow up as in Eq (70). In conclusion, for 
a given laser power, several water vapor lines with a wide range of absorption strengths 
must be selected to accommodate varying atmospheric conditions, i.e. changes in 
humidity and backscatter coefficient, in order to make accurate DIAL measurements of 
atmospheric water vapor over a practical range of ~1- 10 km. 

The next chapter reviews water vapor absorption spectroscopy. Candidate 
absorption lines are selected for DIAL experiments. Linewidth and line center 
wavelength are measured for these candidate absorption lines. Frequency stability 
requirements are determined for the water vapor DIAL laser transmitter and a 
stabilization technique is developed and tested. 


45 



3. H2O VAPOR SPECTROSCOPY, ABSORPTION LINE PROFILING, AND 
LINE LOCKING EXPERIMENTS 

Remote sensing of atmospheric water vapor using the DIAL technique requires 
selecting water vapor absorption lines which have absorption cross sections appropriate 
for the water vapor density and measurement range of interest. For vertical path DIAL 
measurements lines should be insensitive to variations in temperature and pressure 
over the measurement range. Further, criteria must be established for selecting 
absorption lines, taking into consideration not only temperature sensitivity and line 
strength but pressure shift effects, broadening effects, and interference due to 
neighboring lines. Once established, these criteria may be used to select specific water 
vapor lines appropriate for DIAL measurements of tropospheric water vapor. 

The DIAL method employs the difference in the water vapor absorption cross 
section between X on and to estimate water vapor number density in the 
atmospheric path. The absorption cross section Eq.(1) is a function of line strength, 5, 
and linewidth, y, both of which are functions of pressure and temperature. Therefore, 
the temperature and pressure characteristics of S(T) and y(p,T) must be examined 
before selecting absorption lines for DIAL measurements. The three water vapor 
absorption bands in the near infrared, Figure 1, centered at 720, 820, and 940 nm, 
originate in the vibrational ground state with different vibrational overtones as their 
upper states. The weak to moderate line strengths typical of the 720 and 820 nm bands 
are best suited to DIAL measurements in the wetter troposphere while stronger lines 
near 940 nm are better suited for the dryer tropopause and stratosphere. Commercial, 
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high power AIGaAs lasers are most readily available from 780 to 860 nm, therefore the 
820 nm band is best suited to DIAL measurements of tropospheric water vapor using 
AIGaAs lasers. It should be noted that the expressions for water vapor line strength and 
line shape in the 820 nm band are identical to those in the 720 and 940 nm bands. 

Water vapor absorption linewidths are broadened in the atmosphere by two 
processes, pressure broadening and Doppler broadening, which are represented by the 
Lorentz lineshape and the Doppler or Gaussian lineshape, respectively. Variations in 
absorption line shapes due to pressure are not critical since the spectral linewidth of 
AIGaAs lasers (1-10 MHz) is two to three orders of magnitude narrower than typical 
water vapor absorption linewidths (1-5 GHz). Any uncertainties due to pressure shift 
of line center can be minimized by stabilizing the laser wavelength to the center 
frequency of the absorption line at a lower pressure, i.e. higher altitude.^ 5 Therefore, 
pressure shift and broadening effects are not dealt with in this work. 

3. 1 Absorption lineshape, line strength, and cross section 

The shape of a water vapor absorption line can be represented by a Lorentz line 
shape when atmospheric pressure broadening dominates, i.e. altitudes < 1 5 km.^S The 
absorption coefficient for a Lorentz line shape is 


k m (v) = n(J a 



U 2 

+ (r/2) 2 ’ 


(71) 


where n is the number density of water molecules [moled cm 3 ), o a is the absorption 
cross section | cm 2 / molecj, Y is the Lorentz linewidth [cm~^ ), V 0 is the frequency at 
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line center (cm -1 ), and S is the line strength with units of cm~ x j [molecules ■ 
Figure 1 2 shows a plot of the Lorentz line shape using the linewidth and center 
frequency of a water vapor absorption line at 811 .61 7 nm. 



wavenumbers (ctn-1) 

Figure 1 2. Plot of Lorentz lineshape for water vapor absorption line at 81 1 .61 7 nm 
with half-width, y - 0.0837cm -1 / atm. 


The pressure and temperature dependence of the Lorentz linewidth, y(p,T ) , is given 
by46 


Y = 



(72) 


where p is pressure, T is temperature, the subscript "o" denotes initial or STP 
values, and a is the linewidth temperature dependent parameter. Measured values of 
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a, for lines in the 720 nm band, vary widely, from 0.28 to 0.88 with 0.62 the value 
most often used. 47 At line center, the absorption coefficient, k m (v) , reduces to 


K = 


n 
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n -Yo ' 


(73) 


which is the absorption cross section, cr, multiplied by the number density, n. The line 
strength is related to the dipole moment of the water vapor molecule by 4 ® 


. hojL t 
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where the oscillator strength, f , of a dipole transition is 4 ® 
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and v fi is the frequency of the transition, e is the charge of the electron, m e is the 
electron mass, and fi fl is the transition dipole moment. The temperature dependence of 
the absorption line strength, S, for a triatomic molecule, e.g. water, is given by 5 ^ 


S(T) = S(T o y 
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where Q v is the vibrational partition function, Q r is the rotational partition function, 

E" is the ground state energy, h is Planck’s constant, and k is Boltzmann’s constant. 
For the water molecule, the vibrational partition function varies negligibly from unity 
over the temperature range 1 75-3 2 5K and the temperature dependence of the rotational 
partition function is given by 51 


T 

J o 

Kt 
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where j = 3/2 . The expression for the temperature-dependent line strength reduces to 
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Given the temperature dependencies of the line strength and linewidth, the temperature 
sensitivity of the absorption cross section, at line center, may be calculated. 

3.2 Temperature sensitivity of 820 nm water vapor absorption lines 

The temperature sensitivity of water vapor absorption cross section at 820 nm 
is calculated using the same procedure as was previously done for the 720 and 940 nm 
water vapor bands. 46,52 Substituting the temperature dependent expressions for line 
strength Eq (78) and linewidth Eq (72) into the expression for absorption cross 
section, Eq (1 ), yields 
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After simplifying the expression and consolidating terms in the denominator with like 
terms in the numerator, the absorption cross section at line center is written as 
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Taking the derivative of o 0 with respect to temperature once yields the change in 
absorption cross section at line center as a function of temperature, which is written as 
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where the expression within the braces has the form Cj/^fJ-exp^ • t)). Applying 

the product rule for differentiation and simplifying the result yields 
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Recognizing that the expression on the right hand side of Equation (82), excluding the 
curly braces, is equivalent to the expression for absorption cross section, a 0 , as in Eq. 
(80), then the expression simplifies to, 


do 0 _ 
dT 



(83) 


Setting this expression equal to zero yields 
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Since G 0 does not equal zero within an absorption line, the term in the braces on the 
right hand side of Eq. (84) must equal zero for the expression to hold, 

(“~^) + lF = °‘ (85) 

Solving this equation for temperature yields an expression for the temperature neutral 
points of a given absorption line as a function of the ground state energy level, E" , and 
the linewidth temperature dependent parameter, ct, 
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Equation (83) may be plotted for given values of a and E". Dividing da 0 /dT by cr 0 
and plotting the result as a function of temperature yields the fractional change in cross 
section due to changes in temperature. Figures 13, 14, and 1 5 show temperature 
sensitivity of absorption cross section, for Lorentz profile lines, as a function of E" 
for a =0.62, 0.88, & 0.28 respectively. These values of a reflect the most often used, 
0.62, and the measured extremes, 0.88 and 0.28. Therefore Figure 13 may be used in 
general to select temperature insensitive absorption lines while Figures 1 4 and 1 5 may 
be taken as the upper and lower bounds. 



Figure 1 3. The change in absorption cross section, C 0 = sigma, as a function of 
temperature, ground state energy level, E"(cm~ x ), and temperature dependent linewidth 
parameter, a = 0.62. 
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Figure 14 . The change in absorption cross section, o 0 = sigma, as a function of 
temperature, ground state energy level, E"{cm ~ and temperature dependent linewidth 
parameter, a = 0.88. 

The temperature neutral points defined in equation (86) are the zero crossings of each 
E" curve in the figures. For mid-latitude, spring/summer conditions the Earth’s 
troposphere ranges in temperature from 288 K, at the surface, down to 210 K, at the 
tropopause. By selecting water vapor absorption lines with E" = 100 — 250 cm -1 , the 
error in absorption cross section due to temperature variation may be limited to less 
than ±0.1%. In addition to temperature sensitivity, selection criteria must take into 
consideration the absorption line strength and separation from neighboring lines. 
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Figure 1 5. The change in absorption cross section, <7 0 = sigma, as a function of 
temperature, ground state energy level, and temperature dependent linewidth 

parameter, a = 0.28. 

3.3 Absorption line strength and optical thickness 

Optical thickness or optical depth is a concept used to define attenuation of a beam 
propagating through a medium. In our case, the medium is an atmospheric path 
containing water vapor. The path length, x, and volume extinction coefficient, f3 ex , 
which includes scattering and absorption due to both molecules and particles, determines 
the one-way optical thickness, T^, which is defined as 

X 

Tex=\Pex(x)dX- ( 87 ) 
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This is the general expression for the optical thickness, applicable to both homogenous 
paths and paths with variable extinction. A one way optical thickness with a value near 
unity is considered near optimal for making DIAL measurements of water vapor. 53 
Since extinction coefficient and path length determine optical thickness, any practical 
DIAL system will require operation at several different wavelengths in order to access 
absorption lines of varying strength in order to maintain an optical thickness near unity 
over the expected measurement range and water vapor density. The volume extinction 
coefficient (3^ is used interchangeably with k(X,R ) the total extinction coefficient, 

Eq. (38). When the laser is tuned to the center of an absorption line, the total extinction 
coefficient, k(X,R) , is dominated by the water vapor absorption coefficient, k m . 
However, the other attenuation coefficients may not be neglected. For a homogenous 
path, the integral in Eq. (87) reduces to a product, and the optical thickness, as a 
function of absorption line parameters, may be written by substituting Eq. (73) into Eq. 
(87) for yielding 
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X, 
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where n is the density of water molecules along the path of length x. This permits 
calculation of optical thickness for various measurement ranges and selected water 
vapor absorption lines. Using the molar mass of the water molecule, 1 8.01 52 
grams/mole, the mass of a single water molecule is found by dividing the molar mass by 
Avogadro's number, yielding 


Mass(H 7 o) = 18-01 -y r s 3.00 x 10 23 (grams / molecule) . 

2 6.025 xlO 23 


(89) 
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The mean density of water vapor at sea level for mid latitude summer conditions^ j$ 
"1 4.0 g/m3 which corresponds to a water vapor number density of 


n H 2 o = [ h 2°] = — = 4.6 x 10 17 molecules / cm 3 . 

3.00x10 g / molec 


(90) 


Using this number density, the optical thickness for several candidate absorption lines, 
is presented in Table 8 as function of measurement range (horizontal path at sea-level). 


Table 8 


Calculated optical thickness vs measurement range for selected absorption lines. 


Wavenumber Wavelength 

Line Strength 

Linewidth 

Optical Thickness @ Range 

cm-1 

nm 

cm-1 /mol cm-2 

cm-1 *atm 

1 km 

2 km 

5 km 

1 0 km 

12254.533 

816.024 

1.51E-23 

0.0829 

2.70 

5.40 

13.51 

27.02 

12259.498 

815.694 

4.01 E-24 

0.0579 

1.02 

2.05 

5.13 

10.27 

12259.572 

815.689 

1.20E-23 

0.0583 

3.05 

6.10 

15.26 

30.53 

12261.018 

815.593 

6.28E-24 

0.0904 

1.03 

2.06 

5.15 

10.30 

12349.352 

809.759 

2.07E-25 

0.0937 

0.03 

0.07 

0.16 

0.33 

12356.828 

809.269 

4.42E-25 

0.0869 

0.08 

0.15 

0.38 

0.75 

12359.925 

809.066 

1.26E-25 

0.0923 

0.02 

0.04 

0.10 

0.20 

12365.038 

808.732 

4.50E-25 

0.0972 

0.07 

0.14 

0.34 

0.69 

12375.118 

808.073 

1.88E-25 

0.0879 

0.03 

0.06 

0.16 

0.32 


3.4 DIAL absorption line selection criteria 

To reduce errors in the DIAL measurements due to interference from neighboring 
lines, individual absorption lines should be resolved from neighboring lines by 10-20 
pm, or 1-2 absorption linewidths. The following selection criteria have been applied to 
a listing of water vapor absorption lines from the HITRAN database. 

1 ) Single lines separated by 0.2 cm~ x or ~ 10 pm from neighboring lines 
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2) Lower level energies of E" = 100 — 250 cm~ l 

3) Line strengths of 10 -25 - 10 -23 cm -1 / molecule ■ cm~ 2 

Based on the optical thickness calculations in Table 8, temperature sensitivity, and 
wavelength separation, the lines most attractive in the 820 nm band for DIAL 
measurements of atmospheric water vapor are listed in Table 9. 

Table 9 


Water vapor absorption lines which meet selection criteria for DIAL measurements. 


Effective Range 

Wavelength 

(nm) 

Line Strength 
(X 10-23) 

Lower State Energy 
E” (cm* 1 ) 

Short Range (<2 km) 

816.024 

1.51 

300 

815.694 

1.20 

447 


815.593 

0.63 

275 

Long Range (>2 km) 

810.861 

0.017 

142 

81 1.464 

0.045 

173 


808.732 

0.094 

136 


3.5 Absorption line profiling experiments 

An experimental setup has been built using a multipass absorption cell and an 
AIGaAs laser to detect and profile selected water vapor absorption lines near 820 nm. 
This line profile data is used to estimate water vapor absorption linewidth and line 
center wavelength. Figure 1 6 shows a block diagram of the experimental set-up. A 
temperature and current controlled AIGaAs laser is used to profile the absorption lines. 
The laser is used with a chopper wheel, to modulate the light, and lockin amplifier in a 
phase sensitive detection scheme. The wavelength is monitored with a vacuum 
wavemeter, Burleigh WA-20, which yields an absolute accuracy of ±1 part in 10 6 
providing that the laser linewidth is ^ 2 GHz. This corresponds to better than 1 pm 
accuracy at 820 nm but is limited by the instrument readout precision of 1 pm. Once 
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lines have been profiled, the chopper wheel is replaced with a current dither in order to 
perform the derivative spectroscopy measurements. The dither current modulates the 
frequency of the laser, which, when scanned through an absorption line, can be 
demodulated to generate the first, second, or third derivative of the absorption line with 
respect to current. 



Figure 1 6. Block diagram of water vapor absorption line profiling experiment. 


3.5.1 AIGaAs laser diode characterization for line profiling 

The laser diode used to profile absorption lines was a Mitsubishi ML 541 2N 30 
mW device. The laser was fully characterized before it was mounted in the test set up. 
This characterization included measuring the optical output power vs bias current at 
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different case temperatures ( T c ), characterizing the emission spectra with an optical 
spectrum analyzer, and measuring the emission wavelength vs. bias current as a 
function of temperature. Figures 17-20 show PI curves for Mitsubishi device 87- 
1 7003, taken at 7>1 5, 20, 25, and 30°C, respectively. This device, was used to 
probe water vapor lines near 816 nm. The slope efficiency, J] , of the device is estimated 
to be 0.62 mW/mA with a threshold current, l (h , of approximately 42 mA at T c =l 5°C. 

As T c is increased to 20°C, l th increases and r\ decreases to 0.61 mW/mA. The typical 
slope of I th vs T c is 0.66 mA/°C and Tj vs. T c is -0.0014 mW/mA per C 



Figure 17. PI curve for Mitsubishi ML5412N laser diode (serial #87-17003) 
operating at 7’ c = 1 5°C. 
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Figure 1 8 . PI curve for Mitsubishi ML541 2N laser diode (serial #87-17003) 
operating at 7 c =20°C. 



Figure 19. PI curve for Mitsubishi ML541 2N laser diode (serial #87-1 7003) 
operating at T c =2 5°C. 
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Figure 20. PI curve for Mitsubishi ML541 2N laser diode (serial #87-17003) 
operating at r c =30°C. 

The emission spectra of device 87-17003 is shown in Figure 21 for different 
temperatures. These spectra were taken with the laser diode operating at 97 mA bias 
current using a grating optical spectrum analyzer with 0.1 nm resolution. The laser 
operated single longitudinal mode at each temperature. Laser wavelength increased with 
temperature. The -3 dB points on these spectral curves match the instrument 
resolution implying that the laser linewidth is < 0.1 nm. 
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Figure 21. Emission spectra of Mitsubishi ML5412N laser diode (serial #87- 
17003) operating at I op =97.0 mA and T c =15, 20, 25, and 30°C. 

The wavelength of an AIGaAs laser is tunable with bias current. Figure 22 shows 
the current tuning characteristics for different temperatures. This data was taken with 
the grating spectrum analyzer to measure the peak emission wavelength and to ensure 
single mode operation at each temperature. The data for each curve was taken only with 
increasing bias current to avoid hysteresis inherent in AIGaAs laser wavelength tuning. 

The figure also shows the maximum temperature tuning range of this device to be 
approximately 4.0 nm with the laser operating at its 30 mW rated power level. 
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Bias Current (mA) 

Figure 22. Wavelength dependence on bias current for Mitsubishi ML5412N laser 
diode (serial #87-17003) operating at T c = 1 5, 20, 25, and 30°C. 

3.5.2 Profiling water vapor lines 

Profiling water vapor lines and estimating their linewidths and line center 
wavelengths were done by scanning the wavelength of the laser through a known water 
vapor absorption line while detecting a signal transmitted through the White cell with a 
photodiode. The White cell was evacuated and backfilled with pure water to a vapor 
pressure of ~20 Torr. The temperature of the White cell was not controlled but assumed 
constant throughout each measurement. The absorption line profiles were observed by 
modulating the light at ~4 kHz with a chopper wheel and detecting the chopped signal, 
through the White cell, with a photodiode and a Stanford Research Systems SRS 530 
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lockin amplifier. The lockin amplifier acts as a tunable, narrowband filter useful for 
sensing weak signals in the presence of noise. The experimental set-up was controlled 
by a Hewlett-Packard computer which also acted as the data acquisition system. The 
computer recorded the laser diode bias current, lockin error voltage, and wavelength. 
The temperature of the laser diode was sensed initially by an integrated circuit 
temperature sensor (AD-590) which fed a proportional-integral temperature 
controller which controlled the current to a thermoelectric cooler. The temperature 
control loop held the temperature stable to better than 0.1 °C. The AD-590 was later 
replaced by a 1 0 kQ thermistor which improved the temperature sensing resolution to 
0.02°C. 


The laser diode bias current was manually adjusted to its preset starting current 
before each scan. This was necessary to coax the laser into the correct longitudinal mode 
for the desired water vapor line. Once in the right laser mode, the lockin amplifier 
parameters, including sensitivity, dynamic reserve, phase, and time constant, were 
manually set. The bias current was then scanned, under computer control, through the 
water vapor absorption line in 50 or 1 00 fxA steps. At each current step, the computer 

recorded the lockin amplifier error voltage and wavelength. Figure 23 shows the 
absorption signal for a weak water vapor absorption line centered at 81 5.769 nm. 
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Figure 23. Weak absorption feature centered at 81 5.769 nm. 

The vacuum wavelength of this line, recorded in the HITRAN database, is 81 5.772 nm 
and absorption line strength is 2.83 x 10 cm~ x / molecule ■ cm~ 2 . The discrepancy in 
observed wavelength may be attributed to either a vacuum error in the wavemeter, 
which would tend to shorten the wavelength due to an increase in refractive index inside 
the interferometer, or the pressure shift effect of line center due to the partial pressure 
of water vapor in the White cell. This data was taken with an absorption cell path length 
of ~ 20 m, vapor pressure < 20 Torr, and laser diode temperature of 1 9°C. The 
expected absorption for this line strength, vapor density, and path length is <0.01%. 
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A stronger absorption line was scanned and is shown in Figure 24. The vacuum 
wavelength of this line is 81 6.024 nm. The observed wavelength however was 81 6.027 
nm. The absorption strength for this line is 1.51 x 10 ~ 23 cm- 1 /molecule cm~ 2 . The 


expected absorption for this line 


strength, vapor density, and path length is ~5%. 


0.40 


816.035 




Current (mA) 

Figure 24. Profile of a strong absorption line centered at 81 6.024 nm. 

The wavelength is also plotted in this figure. A simple line fit has been performed on the 
wavelength and its equation included. The slope of the fit indicates a wavelength tuning 
rate of approximately 8 pm/mA. 

The linewidth and line strength of an absorption line may be estimated from 
measurements made using second harmonic spectroscopy techniques. 55, 56 These 


67 




measurements used a low amplitude dither current to modulate the frequency of the laser 
emission, a silicon photodiode to detect the light transmitted through an absorption cell, 
and a lockin amplifier to observe the first, second, and third derivatives of the 
absorption signal. 

The chopper wheel was removed and a 4 kHz, 2 mVpp sine wave was superimposed 
onto the laser via a 47 Cl series resistor. The 47 Cl resistor in series with the laser 
diode creates a ~50 Cl load. The 2 mV sine wave translates to a ~40 jiA dither current 
which modulates the center frequency of the laser. The laser beam was directed through 
the White cell and detected by an EG&G Si PIN photodiode. The photodiode signal was 
demodulated by the lockin amplifier with an integer multiple of the dither frequency, /. 

It has been shown that demodulating with a integer multiple off yields the integer 
derivative, i.e. demodulating with / yields the first derivative, demodulating with 2/ 
yields the second derivative, etc. Figure 25 shows the first derivative of the 81 5.769 
absorption line shown in Figure 23. 
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Figure 25. First derivative of water vapor absorption line at 81 5.769 nm. The 
wavelength is also shown. 

The signal strength for the first derivative curve in Figure 25 is higher than the 
absorption signal in Figure 23 because neutral density filters were removed from the 
beam. These neutral density filter were used to attenuate the beam prior to entering the 
White cell so as to not saturate the photodiode. The first derivative of the 81 6.024 nm 
absorption line is shown in Figure 26. The wavelength has also been plotted in this 
figure. 


69 


Wavelength (nm) 





816.032 

816.030 

816.028 

816.026 

816.024 

816.022 

816.020 

816.018 


a> 

< 

<D_ 

CD 

3 

CO 


3 

3 


Current (mA) 


Figure 26. First derivative of water vapor absorption line at 81 6.024 nm. Also 
shown is the wavelength vs. bias current. 


Demodulating the absorption signal with 2/ yields the second derivative of the 
absorption signal with respect to current. Figures 27 and 28 show the second derivative 
of the absorption lines centered at 81 5.769 and 81 6.024 nm, respectively. The second 
derivative curves are all inverted in amplitude. The center of the second derivative 
should be a positive valued peak reflecting the steep positive slope of the first 
derivative. This problem has been traced to the autophase (AP) command sent to the 
lockin amplifier during the initialization and subsequently corrected. Estimation of 
linewidth and wavelength are unaffected by the relative phasing of the lockin however. 
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Figure 27. Second derivative with respect to current of the absorption feature 
centered at 81 5.769 nm. 

The full width at half maximum (FWHM) of an absorption line may be estimated 
from the second derivative curve. The zero crossings of the central peak of the second 
derivative correspond to the peaks of the first derivative curve which in turn 
correspond to the points of steepest slope (inflection points) on the absorption line. 
These points are close to the half maximum points, as long as the absorption is small, < 
5%, and therefore can be used to define the linewidth. From Figures 27 and 28, the 
current linewidths are similar and estimated to be approximately 0.350 mA. Using the 
experimentally determined current tuning ratio of 8 pm/mA the lines are 
approximately (0.350mA)*(8pm/mA) = 2.8/?m wide or lApm halfwidth. Converting 
the halfwidth from pm to cm~ x yields a measured halfwidth of 0.021 cm~ x . At 20 Torr 


71 




total pressure, the lineshape is dominated by Doppler broadening and the Doppler 
halfwidth of this line may be calculated using 


«D = 


[ 2kln2 if 


( 91 ) 


where v 0 (cm ~ x ) denotes line center, k(ergs / K) is Boltzmann’s constant, T(K ) is the 
temperature, c(cm / sec) is the speed of light, and m(g / molecule) is the mass of the 
water molecule. For the line at 81 6.024 nm the Doppler halfwidth is 


a D — 


12,254cm- 1 


3 x 10 10 cm / sec 


2(1.38 xl0 16 cr^//ir)ln 2- 


29 A K 


1>2 


3 x 10 ~ 23 g / molecule J 


(92) 


or a D = 0.0177cm" 1 . The difference between measured and calculated halfwidths is 
1 8.8%. This difference is attributable to the choice of the zero crossings of the second 
derivative curve to define the linewidth. The peak absorption should be estimated from 
the absorption line profile, Figures 23 & 24, and the half-absorption points 
determined. However, for weak absorption, <5%, the zero crossings of the second 
derivative curve are an acceptable approximation to the FWHM. 
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Figure 28. Second derivative with respect to current of the absorption feature 
centered at 81 6.024 nm. 


Demodulating the absorption signal with 3/ yields the third derivative of the 
absorption feature with respect to current. Figures 29 and 30 show third derivatives of 
the absorption lines centered at 81 5.770 and 81 6.024 nm. The third derivative is 
desirable for use as the frequency discriminant in the line locking electronics due to its 
long linear region through the line center and its nearly zero amplitude offset. The first 
derivative includes a nearly constant amplitude offset due to the positive slope of the 
absorption signal which is due to the linear intensity increase with laser diode bias 
current. This amplitude offset moves the zero crossing of the first derivative curve 
away from the absorption line center. Therefore locking the laser to the zero crossing of 
the first derivative would cause an underestimation of water vapor density, in an 
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Figure 30. Third derivative of the absorption feature centered at 81 6.024 nm. 

The requirements for closed loop frequency locking electronics have been 
established from these water vapor absorption line profiling experiments. Pressure 
broadened lines have Lorentzian lineshape Eq. (71) which determines how the 
absorption coefficient, k m , changes as the laser frequency varies from line center. To 
ensure a change of less than 1 0% in the absorption coefficient, the laser must be locked 
to 0.34 of the Lorentz halfwidth. For the Mitsubishi lasers, with ~8pm/mA wavelength 
tuning rate, this corresponds to a bias and feedback current stability of ~0.2 mA. 

The 1 00 mW AIGaAs lasers procured for the AIGaAs lidar transmitter (Spectra 
Diode Labs SDL 5410) were specified at 81 6±5 nm. One laser at a time was mounted in 
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a laser diode header and collimated. Each laser was then checked with a current source 
and a fiber coupled optical spectrum analyzer. Each laser operated in a single 
longitudinal mode up to 100 mW and was centered in wavelength at ~81 2 nm. This 
required identifying several new candidate absorption lines at 812 nm and profiling 
these lines with the same water vapor spectroscopy set-up, Figure 1 6. Two moderately 
strong lines in this region, 811 .006 and 81 1 .61 7 nm, were selected for use later in the 
horizontal path water vapor DIAL measurements. These lines permitted using a shorter 
absorption cell path length. Consequently, a 30 cm long cell with angled windows was 
used in a multipass configuration. The cell was evacuated and backfilled with ~ 20 Torr 
of water vapor. The temperature of the cell was not controlled. 

A first derivative scan of a line at 81 1.617 nm is shown in Figure 31. The line 
strength of this absorption line is 2.54 x 10 -24 cm -1 / molecule • cm~ 2 and it has a 
pressure broadened linewidth of 0.0837cm _1 . Also shown in the figure is the wavelength 
vs. bias current curve. The SDL-5410 laser exhibits a slower current tuning rate, 

3.3pm / mA, than the Mitsubishi ML5412N laser, 8 pm/mA. This line was 
subsequently used in both the horizontal integrated path and the range-resolved water 
vapor DIAL measurements. 
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Figure 31 . First derivative scan of 81 1 .61 7 nm absorption line using 90 cm 
absorption path length in a 30 cm absorption cell. 

The current linewidth measured from Figure 3 1 is 0.7 mA which corresponds to a 
linewidth of (0.7mA) ■ (3.3pm / mA) = 2.3 pm and a halfwidth of 0.0175 cm -1 . This 
compares well with the calculated Doppler halfwidth of a D = 0. 0177cm -1 . The slope of 
the central portion of the of the first derivative curve is 

-^^ = 4.28 mV/mA. (93) 

0.1mA V 1 
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3.6 Locking an AIGaAs laser to a water vapor absorption line 


An AIGaAs laser has been locked, in frequency, to an absorption line using a 
computer to close the feedback loop. Using a setup similar to the line profiling 
experiments, the AIGaAs laser was dithered with a 3 kHz, 0.4-0.8 mA sine wave. The 
collimated laser beam was transmitted through an absorption cell and a photodiode 
detected the absorption signal. The lockin amplifier demodulated the absorption signal at 
3kHz. This generates the first derivative of the absorption signal with respect to 
current. The wavelength of the laser diode is controlled by a feedback current which is 
summed, at the laser, with the bias current. The feedback current is determined by the 
lockin voltage (which is manually set to zero at line center) and by the loop gain which 
is set by a computer program controlling the feedback current source. The software, 
which attempts to zero the lockin error, adjusts the gain by examining the lockin error 
voltages for the last five points. Figure 32 shows a flow chart of the computer algorithm 
which was used to stabilize the laser wavelength to the water vapor absorption line. The 
absolute value of the sum of five samples of the lockin voltage is compared with the sum 
of the absolute values of the same lockin voltages. If the absolute value of the sum is > 
the sum of the absolute values, the loop gain is increased. Conversely if the absolute 
value of the sum is less than the sum of the absolute values, the loop gain is decreased. 
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Figure 32. Algorithm for laser stabilization to water vapor absorption line. 


Frequency stability of the stabilized laser has been measured and recorded using 
the White cell and lockin amplifier. By ensuring that the laser frequency is stable over 
the typical lidar measurement period (1-5 minutes), DIAL measurement errors due to 
wavelength drift of the laser will be reduced. The error signal, generated by the lockin 
amplifier sensing the transmission signal through the absorption cell, is proportional to 
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the fluctuations of the laser frequency. 58 Figure 33 shows the lockin error signal as a 
function of time while the laser was actively stabilized to the center of the absorption 
line at 81 1 .61 7 nm. Each data point corresponds to approximately 1 second. The lockin 
error voltage fluctuates approximately +20 / -30/J.V over the ~200 second measure- 
ment period. From Figure 31 the amplitude (peak to peak voltage) of the first 
derivative of the absorption line at 81 1 .61 7 nm is measured to be 3 mV. The slope of 
the central portion of this curve is estimated to be 4.28 mV / mA. 



data point 


Figure 33. Scan of lockin voltage and laser wavelength while laser was actively 
stabilized using the computer feedback loop. Each data point corresponds to ~1.0 second. 
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The fraction of the absorption linewidth over which the laser fluctuated during the 
measurement in Figure 33 is the ratio of the lockin error fluctuations to the amplitude 
of the first derivative, 


50pV 

3mV 


0.016, 


or one part in 60. The center frequency stability of the laser is 


0.05/nV „„ , a a 

Y-p ... - ' 3. 3 pm / mA = 0.04pm , 
4.28 mV / mA 


(94) 


(95) 


over ~200 seconds. 

The demonstrated frequency stability of the actively stabilized laser was 
sufficient to attempt DIAL measurements of atmospheric water vapor. The next chapter 
describes the hardware assembled for two prototype AIGaAs lidar systems, a single color 
cloud and aerosol lidar and a water vapor DIAL system. 
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4. PROTOTYPE AIGaAs LIDAR SYSTEMS 


Two prototype AIGaAs lidar systems have been assembled and tested. The first is 
a single color cloud and aerosol lidar and the second is a water vapor DIAL system. Each 
system is based on a single 1 00 mW AIGaAs laser diode, a 20 cm diameter Schmidt- 
Cassegrain telescope, and a silicon Geiger-mode avalanche photodiode {APD). The 
primary difference between the two systems is the method used to modulate the intensity 
of the laser diode emission. Both systems use the same receivers and detectors as well as 
the same histogrammer, data acquisition software, and correlation software. The AIGaAs 
laser diode used for both lidar systems is discussed in this chapter followed by a detailed 
description and diagram of each transmitter subsystem. A complete description of the 
receiver common to both lidar systems is presented as is a description of the data 
acquisition and histogramming hardware and software. 

4. 1 AIGaAs laser diodes 

The laser used in both transmitters is a Spectra Diode Labs SDL-5410-G1 
AIGaAs diode laser. This is a single element index guided AIGaAs laser and has a nominal 
linewidth, inversely proportional to optical power, of 10-100 MHz. Linewidths of 
AIGaAs lasers have been narrowed to 330 kilohertz with electrical feedback 59 . Single 
element laser diodes tend to operate, spectrally, in a single longitudinal mode, although a 
small percentage of power is evident in neighboring modes. The amplitude of these side 
modes is typically 20-30 dB lower than the dominant mode. A plot of output power vs. 
bias current is shown in Figure 34. This PI curve shows the optical power for a laser 
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diode bias current between 50 and 1 00 mA. The slope of the line is very nearly 1 W/A 
which is referred to as slope efficiency. The laser diode package is shown in Figure 35. 



Figure 34. PI curve for SDL 5410 G1 AIGaAs laser diode. 

The emission wavelength of semiconductor lasers is tunable in both temperature 
and bias current. As mentioned in Chapter 3, typical tuning coefficients are 70 pm/°C 
and 3.3 pm/mA. These tuning characteristics are used to stabilize the lidar transmitter 
to the center of an absorption line for the water vapor DIAL system. The emission 
aperture of a single element laser diode is ~ 1x3 fim and light diverges rapidly with a 
3:1 aspect ratio into a cone ~10°x30°. This requires using collimating optics with a 
high numerical aperture to collect as much light as possible into a usable beam. Two 
commercially available collimating lenses were used. A 0.5 numerical aperture lens 
was used in the single color cloud and aerosol lidar. This lens, the LDCO-53-N from 
Optics Plus Inc., has three spherical elements with special multilayer anti-reflective 
(AR) coatings on all air-glass interfaces. A 0.55 NA glass asphere from Corning, model 
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3501 50, was used in the water vapor DIAL system. This singlet had broadband AR 
coatings on both sides. 



On 

L 


8j. 


Figure 35. Laser diode in the SOT-148 window package [dimensions in inches(mm)] 

The AIGaAs lasers used are single element, index guided AIGaAs lasers. The 
Mitsubishi ML541 2N used in the water vapor spectroscopy experiments are 30 mW, 
single mode lasers and were selected for their 81 5-822 nm emission wavelength and 
their wavelength stability. The Spectra Diode Labs SDL 5410 used in the lidar systems 
are 1 00 mW, single mode lasers and were selected for their high average power, 
ruggedness, and reliability. 

When the wavelength stability of the laser is not crucial, i.e. for aerosol and 
cloud lidar measurements, direct current modulation of the laser diode is used to 
modulate the intensity of the outgoing beam. However, for DIAL measurements, the 
frequency deviation of the laser wavelength caused by modulating the drive current 
precludes maintaining the frequency stability required for these measurements. For 
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this reason, an external, electro-optic modulator was used to modulate the intensity of a 
frequency stabilized CW AIGaAs laser. 

4.2 Single color cloud and aerosol lidar transmitter 


In the single color cloud and aerosol lidar, the current to the laser diode is 
modulated with a maximal length PN code. This current modulation directly modulates 
the intensity of the laser output. There is no monitoring of the laser diode wavelength. 
A system diagram is shown in Figure 36 and system components are listed in Table 10. 
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Figure 36. Cloud and aerosol AIGaAs lidar system diagram. 


The laser transmitter for the single color cloud and aerosol AIGaAs lidar is based 


on the 1 00 mW laser diode, Spectra Diode Labs SDL-541 0-G1 . The cap which protects 


the laser diode has been removed so that a short focal length collimating lens can be used. 
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An optical and mechanical mount for the laser diode and collimating lens was designed 
previously 6 ® and is used in the AIGaAs lidar transmitter. The laser diode header 
provides a rugged optomechanical design, good thermal interface, and an electrical 
interface capable of > 100 MHz modulation. 


Table 10 

Cloud and aerosol lidar system components. 


AIGaAs laser 
Collimating lens 
Beam expander 
PN code generator 
Receiver telescope 
Optical bandpass filter 
Detector 
Discriminator 
Histogrammer 
System computer 


Spectra Diode Labs, SDL 5410-G1 
Optics-Plus LDCO-53-N, NA=0.5 
Rodenstock model #4401 2x-8x, adjust. 
Custom, 1 MHz, 31-4095 bit codes 
Meade 20 cm diameter F/6.3 
Andover, 810 nm, 10 nm FWHM 
EG&G SPCM-200 Photon Counting Module 
Tennelec, TC-453 

Custom, 4095 bits long, 32 bit depth 
Gateway 486/33 


4.2.1 Current and temperature controller 


The single color cloud and aerosol lidar requires only ±0.1 mA bias current 
accuracy and ±0.1 °C temperature stability, therefore an integrated current- 
temperature controller was used in conjunction with a digital modulator circuit to 
control the laser diode. The proportional-integral temperature controller works with 
an Analog-Devices AD-590 integrated temperature sensor and a thermo-electric cooler 
(TEC) to form a closed loop temperature controller with 0.1 °C accuracy. The current 
controller is manually adjusted to a bias current up to 60 mA with 0.1 mA accuracy and 
a modulation current up to 1 00 mA. The digital modulator circuit accepts an ECL digital 
signal and converts it to an analog modulation current which is summed with the bias 
current at the diode. 
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4. 3 Water vapor DIAL transmitter 


The water vapor DIAL system, using the same SDL-5410-G1 laser diode, 
requires that the laser wavelength be tuned in wavelength both to the center of a water 
vapor absorption line and to a nearby non-absorbing region. To achieve this tuning and 
the necessary current stability, two programmable current sources are used. A 
Hewlett-Packard computer was used to control both current sources. A separate 
temperature controller was used in a closed loop configuration to stabilize the tem- 
perature of the laser diode header to within 0.02°C. An electro-optic light intensity 
modulator was used to impress the PN code onto the laser intensity. Therefore, the laser 
is run CW and no digital modulation circuitry is needed. A smaller laser diode header 
was designed for the SDL-5410-G1 laser package, Figure 35, and single element 
aspheric collimating lens. The laser diode cap including window, which covers the laser 
diode, was removed so that the Corning asphere lens could be used. A system diagram is 
shown in Figure 37 and major system components listed in Table 1 1 . 
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Figure 37. AIGaAs lidar water vapor DIAL system diagram. 

The transmitter for the water vapor DIAL system has more stringent wavelength 
requirements than the cloud and aerosol lidar. Therefore, the integrated current and 
temperature controller previously used for the cloud and aerosol lidar was replaced by 
two programmable current sources (Keithley model 220 programmable current 
source) running in parallel and a proportional-integral temperature controller 
(Alpha-Omega TC Series 2 temperature controller). 
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Table 1 1 


Water vapor DIAL system components. 


AIGaAs laser 

Laser current sources 

Beam expander 

Collimating lens 

Xmit temperature controller 

Xmit temperature sensor 

Active temp, device 

Absorption cell 

Lockin amplifier 

PN code generator 

External light modulator 

Receiver telescope 

Optical bandpass filter 

Detector 

Discriminator 

Histog rammer 

Data acquisition computer 


Spectra Diode Labs, SDL 5410-61 
Keithley 220 program, current source 
Rodenstock model #4401 2x-8x, adjust. 
Corning molded glass asphere, NA=0.55 
Alpha-Omega TC Series 2 
1 0 kfl thermistor 

Marlow, SP-1666 Thermo-electric cooler 
Connecticut Instruments, Pathfinder Cell 
SRS Model 530 Dual phase lockin amp. 
Custom, 1 MHz, 31 - 4095 bit codes 
ConOptics Model 3, Electro-Optic Modulator 
Meade 20 cm diameter F/6.3 
Andover, 810 nm, 5 nm FWHM 
EG&G SPCM-200 Photon Counting Module 
Tennelec, TC-453 

Custom, 31-4095 bit length, 32 bit depth 
Gateway 486/33 


4.3.2 Current sources and laser frequency control loop 

The two programmable current sources were controlled by a Hewlett-Packard 
9000 series computer via the IEEE-488 computer interface. One current source was 
used to bias the AIGaAs laser to its CW operating point while the second current source 
was used to provide the closed-loop feedback current to the laser. This feedback current 
was set by a computer program which senses the error signal, approximately once a 
second, from the lockin amplifier The wavelength of the laser diode is controlled with 
the feedback current using the 3.3 pm/mA transfer function of the laser previously 
measured in Chapter 3. The feedback current is determined by computer algorithm, 
Figure 32. 
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4.3.3 Temperature controller 


The temperature control loop uses an proportional-integral temperature 
controller, a 1 0 k£2 thermistor, and a custom thermo-electric cooler (TEC) to close the 
loop. The use of a 1 0 k£2 thermistor with the temperature controller provides 0.02°C 
resolution in sensing temperature. The temperature controller is bipolar, capable of 
driving as well as sinking current. The TEC device, which uses the Peltier effect, is also 
a bipolar device which can heat as well as cool. The Peltier cooler uses many junctions 
of two dissimilar metals, which generates a heat flow at the junction when a voltage 
potential is placed across it. 

4.3.4 Electro-optic light modulator 

The EO modulator, Con-Optics Model 3 Light Intensity Modulator, employs a 
Potassium Dideuterium Phosphate (KD*P) electro-optic crystal to digitally modulate 
the polarization vector of the laser beam and a Glan-Thompson polarizer, oriented at 
90" to the input polarization, to convert the polarization modulation into intensity 
modulation. This technique is attractive for AIGaAs lasers since they have polarization 
ratios greater than 1 00: 1 . The EO crystal is electrically biased to the half-wave voltage 
and the digital TTL modulation signal rotates the polarization between maximum 
extinction and maximum transmission. A linear polarizer or analyzer on the input 
improves the performance of the modulator by initially rejecting any unpolarized light. 
KD*P was chosen because of its > 90% transmission at 820 nm and its 80% modulation 
depth at a modest (<100VDC) half-wave voltage. 

The modulator has been characterized with the transmitter laser using a fast Si 
APD detector and power meter. The optical power and modulation depth were first 
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optimized by adjusting alignment, rotation, and bias voltage of the modulator. The power 
was measured at the output of the laser diode header, after the beam splitter, at the input 
of the EO modulator, and again at the output of the EO modulator. This data is shown in 
Figure 38. Optical PN code waveforms were recorded at the output of the EO modulator 
and are shown in Figure 39. The input power to the modulator was 5. 1 7 mW and the 
average output power from the modulator was 1.71 mW. The peak transmission is 
therefore 66% assuming a 50% duty cycle of the PN code modulation. 
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Figure 38. Power vs current characterization of the AIGaAs DIAL transmitter. 
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Figure 39. PN code optical waveform produced by the EO modulator with a -1 25 VDC 
bias. 


With the laser blocked, the detector output was -1.125 mV. The digital zero and one 
light levels were -2.55 mV and -9.875 mV respectively. The modulation depth is 


(-9.875ffiV)-(-2.550mV)| ^ 

8. 750m V 


(96) 


for a modulation depth of 84%. 
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4.4 PN code generator 


In preparation for this work, shift register code generator circuit designs were 
reviewed and a general purpose maximal length sequence generator was designed. The 
design was implemented in TTL fast logic (74F-series) to ensure good quality modula- 
tion waveforms up to 10 MHz bit rates. A schematic diagram of the PN code generator is 
shown in Figure 40. Twelve D flip-flops form the shift register elements of the PN code 
generator. Selected outputs are fedback to the input of the first shift register element 
via a network of exclusive-OR gates. 

The code generator is capable of driving ECL, TTL, and 50 Q loads and provides a 
synchronization pulse at the start of every code sequence. The code generator produces 
sequence lengths from 31 bits to 4095 bits. To test that the circuit correctly produced 
maximal length PN codes, a sequence at each code length was accumulated and cross- 
correlated with an ideal code generated by a computer algorithm. Each code sequence 
produced by the code generator correlated perfectly with the computer generated codes. 
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Figure 40. Schematic diagram of PN code generator with adjustable code lengths from 
31 to 4095 bits. 


4.5 Beam expander & pointing mirror 

The laser diode beam divergence must be known and less than the receiver FOV. 
The diffraction-limited divergence, 6 , of a TEM 0 o mode laser, measured at the 1 / e 2 

points, is 
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where X is the wavelength and d Q is the diameter of the beam waist. An adjustable beam 
expander is used to increase the outgoing laser beam waist, d a , in order to reduce the 
laser divergence to less than 100 prad. The beam expander, Rodenstock model #4401 
2x-8x expander, has a maximum exit aperture of 31 mm at 8x expansion. This 
corresponds to a diffraction limited divergence of ~3.5 prad. A precision mirror mount 
is used to direct the outgoing beam to within the receiver FOV. 


4.6 Receiver 


The receiver system used for both the single color cloud and aerosol lidar and the 
water vapor DIAL system is based on a 20 cm telescope, additional aft optics, bandpass 
filter, and a silicon avalanche photodiode used in the photon counting or Geiger mode. 

4.6.1 T elescope/receiver optics 

The lidar receiver uses a commercially available, Meade 20 cm diameter, F/6.3 
Schmidt-Cassegrain telescope. A 200 pm pinhole was placed at the primary focus of the 
telescope creating a field stop. The receiver field of view (FOV) is 


FOV = 


200/im 
1280 mm 


~ 156j urad. 


(98) 


Light passing through the field stop is recollimated by a 50 mm focal length achromat. 
The field stop and achromat are mounted to the rear of the telescope in an adjustable 
black delrin housing. Once assembled, the telescope, including recollimation optics, was 


mounted to the azimuth/elevation (AZ/EL) stage of a 40 cm white light collimator. The 
collimator filled the aperture of the receiver telescope with plane waves of white light. 

The pointing of the telescope was adjusted until light passed through the field stop. The 
focus of the telescope was adjusted to maximize light through the field stop. Telescope 
pointing was checked and adjusted often during this process. Once the telescope focus was 
optimized, the aft optics were adjusted to provide collimated light to the detector module. 

A small autocollimating telescope was used to monitor the output of the lidar telescope 
and aft optics to ensure that collimation was achieved. Once collimated, the aft optics 
were marked and the threads of the adjustable delrin housing tacked in place with 
adhesive. 

4.6.2 Si Geiger-mode APD/photon counting detector 

The detector used was a silicon avalanche photodiode (Si APD) biased in the 
Geiger-mode of operation. ^ In this mode, single photons may be detected as electrical 
pulses. An incoming photon is absorbed and produces a photoelectron which is amplified 
in the avalanche process. A circuit compares the pulse voltage with a preset level and 
outputs a TTL signal when the pulse voltage exceeds the preset level. One drawback of the 
Geiger-mode APD is the dead time between output pulses. The device, which is passively 
quenched, requires approximately 200 nsec to recharge. Hence, it is limited to count 
rates below 10 6 counts / sec for a linear response and 10 7 counts / sec using a correction 

factor. 

The maximum sensitivity of this detector is confined to a radius of 50 pm from 
its center which makes the receiver optical design critical. A 20 mm focal length laser 
doublet focuses the light onto the active area of the detector. To reduce the background 
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photoelectron count rate, an optical bandpass filter was mounted in front of the focusing 
doublet. The bandpass filter and laser doublet fit in an adjustable aluminum housing 
which mounts to the front of the photon counting module. This black anodized housing 
effectively seals the detector from stray light. The assembled detector module was 
mounted to the AZ/EL stage of the 40 cm white light collimator and adjusted while 
monitoring the photoelectron count rate on a frequency counter. The pointing of this 
module was sensitive due to the small active area of the detector and the short focal 
length lens. Therefore, pointing was checked and adjusted whenever the focus of the 
detector module was adjusted. Once optimized, a knurled ring was used to lock the focus 
adjustment. 


4.6.3 Alignment of integrated receiver 

The telescope, with its integral aft optics, and detector subsystem were 
integrated to a common aluminum breadboard. The optical alignment of these two 
subsystems is critical. The telescope and field stop define the lidar receiver optical axis 
and the detector subsystem must be aligned to it. The integrated receiver was mounted to 
the AZ/EL stage of the 40 cm white light collimator. Initially, the detector was blocked 
with a small piece of cardboard to prevent saturation. The AZ/EL of the collimator was 
adjusted to maximize light through the receiver optics. The detector position was 
adjusted visually to the center of the recollimated light. The collimator light intensity 
was then attenuated, the detector turned on, and the photoelectron count rate monitored 
on a frequency counter. Detector position was adjusted to optimize the photoelectron 
count rate. Once the detector was bolted in place no further adjustments to the receiver 
were necessary. 
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4.7 Data acquisition & histogramming 


A histogramming circuit, shown in Figure 41 , was developed ^7 which 
accumulates received photoelectrons into memory elements synchronously with the 
transmitted PN code sequence. There is one memory element per bit of the PN code 
sequence and the unit is capable of accumulating counts at a 2 MHz count rate 
continuously with no dead time. The histogrammer was built as a single AT expansion 
board to fit within the Gateway 486 system computer. Software has been developed to 
control the histogrammer, transfer data from the histogramming memory to the 
computer hard drive, and compute the cross-correlation function. 
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4.7.1 Histogrammer operation 


The histogrammer and PN code generator are clocked by an external frequency 
generator. The histogrammer senses an 8 bit long sync pulse from the PN code generator 
which determines the start of the histogram. The sync pulse resets the address counter 
in the histogrammer. The address counter begins counting and continues counting until 
the next sync pulse. Meanwhile, a pulse catching circuit consisting of several D flip- 
flops latches an in-coming photoelectron pulse generated by the detector and 
simultaneously latches the address, where the detected pulse occurred, into an address 
buffer. The memory contents at this address are moved to a photon adder where its value 
is incremented by one. The new value is moved back to the data buffer and written to the 
same memory address. Two parallel memory units, each 4 Kbits long by 32 bits deep, 
are used to eliminate any dead time associated with transferring the memory contents to 
the computer. While one memory unit is accumulating photoelectron counts the other is 
writing data to the computer. When the counting memory unit gets to half full, a flag is 
set and the computer disables this unit and enables the empty memory unit to start 
counting incoming pulses. The computer then begins transferring the memory contents 
of the nearly full memory unit to the hard drive. 

Software has been developed to control the histogramming electronics as well as 
manage the accumulated histograms, compute correlation functions, and display the 
histograms and correlations on the computer display. 

The histogrammer was tested by applying clock and sync signals while holding the 
input to the pulse catching circuit first at ground and then to a high TTL signal. The 
histogrammer was allowed to accumulate for several seconds. With a grounded input, the 
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histogrammer performed as expected and recorded no counts in any histogrammer bin. 
With the input tied high however, the histogrammer did not perform perfectly. There 
was typically a discrepancy of one count per bin between the bins in the first half of the 
histogram and the bins in the second half. Also, the last bin of the histogram would not 
accumulate counts. The one count per bin discrepancy had a negligible effect on 
correlation functions of long integration time measurements and was ignored. The zero 
count in the last bin of the histogram was found to have significant impact on the 
correlation function, especially during long integration periods or during measurements 
with high average count rates. For this reason the system software was modified to set 
the count in the last bin to the average count in all of the histogrammer bins. While this 
is not correct, simulations have shown that it does not induce a significant error in the 
correlation function and that it does reduce the error caused by the zero count in the last 
bin. 


The next chapter presents lidar measurements made with the two prototype 
AIGaAs lidar systems. The lidar data is analyzed and compared with theory wherever 
possible. 
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5. LIDAR MEASUREMENTS & DATA ANALYSIS 


Once assembled and boresighted, the AIGaAs lidar was used to make several 
different lidar measurements. These lidar measurements consisted of ranging over 
horizontal paths to terrestrial targets and measurements of clouds and aerosols over 
slant paths. In addition, DIAL measurements of atmospheric water vapor were made over 
horizontal paths in both integrated path and range resolved modes. This chapter 
summarizes results from these measurements and provides analysis of the data, 
including range correction, SNR calculation, and an estimate of laser diode frequency 
stability. Also, an error analysis of the water vapor DIAL measurements is performed. 

5.1 Horizontal path measurements to terrestrial targets 

The AIGaAs lidar was used to range over horizontal paths to terrestrial targets at 
ranges of 0.5 to 13 km. The lidar was operated at a single, unstabilized wavelength with 

direct current modulation by a PN code. The AIGaAs lidar is a biaxial lidar system with a 
small transmitter divergence, < 100 /urad, and a narrow receiver FOV, 160 firad, it is 

therefore susceptible to small changes in the boresight alignment (the transmitter to 
receiver alignment). The overlap function was optimized for ranges >5 km. This was 
done by adjusting the transmitter pointing mirror while observing photoelectron count 
rates from distant targets. The strong backscattered signals from these targets reduced 
integration times and permitted near real-time correlation of the received signal with 
the transmitted signal which facilitated the alignment process. By comparing the 
backscattered photoelectron count rate with the expected count rate, calculated using 
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the single scattering lidar equation and the system parameters, any serious problems 
could be identified. 


5.1.1 Alignment method #1 

Two methods were used to align the transmitter and receiver axes of the AIGaAs 
lidar. The first boresight method used a target board on which cross-hair marks were 
made at positions with the same vertical and horizontal separation distances as the 
transmitter pointing mirror and receiver telescope. A 20 mW HeNe laser was mounted 
to the lidar transmitter to provide a visible (632.8 nm) laser pointer. A third mark 
was added to the boresight board corresponding to the HeNe’s position. The telescope was 
focused to its respective mark on the boresight board by viewing through the receiver 
telescope with either an eyepiece or a charge coupled device (CCD) camera. The HeNe 
laser and AIGaAs laser transmitter were turned on and their pointing adjusted until each 
laser overlapped its proper mark on the boresight board. Using a CCD camera made the 
AIGaAs laser spot visible on the video display. The boresight board was then moved 
farther away from the lidar while the position of the boresight marks and laser spot 
were monitored. The telescope was refocused and the AIGaAs laser transmitter pointing 
mirror adjusted to bring it back into alignment. Finally, the telescope was focused at 
objects on the local horizon, 12-15 km away. This insured that the telescope was 
focused effectively at infinity. 


5.1.2 Alignment method #2 

The second alignment method employed a 1 2.7 cm hollow corner cube retro- 
reflector (CCR) to redirect the transmitter laser light back into the telescope. With the 
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AIGaAs laser operating above laser threshold, the transmitted light was directed into the 
edge of the CCR. The manufacturing tolerance on the CCR was 2 arcsec (9.7 jirad). 

Since the CCR translates the return beam by approximately its aperture size it allowed 
the transmitter to reflect directly into the telescope. A piece of translucent mylar was 
placed at the exit aperture of the telescope. With an infrared (IR) viewer, the 
transmitter pointing mirror was adjusted until a bright spot appeared on the mylar. 

The intensity of this spot was maximized with further fine adjustment of the pointing 
mirror. 


With adequate neutral density filters in the transmitter path, the photon- 
counting detector was turned on and the photoelectron count rate monitored with a 
frequency counter. The photon counting detector had been previously aligned to the axis 
of the telescope using the 1 6” white light collimator. The transmitter pointing mirror 
alignment was adjusted to maximize the count rate. Once the alignment was optimized, 
the laser drive current was reduced to less than 1.0 mA and the neutral density filters 
were removed. This was to check for angular misalignment due to using multiple neutral 
density filters. Since the optical system alignment was preserved, the need for ND 
filters was eliminated and in all subsequent experiments this alignment method was used 
for boresighting the AIGaAs lidar. 

5.1.3 Measurements to terrestrial targets 

Several horizontal path lidar measurements to terrestrial targets were made 
from the rooftop lidar laboratory of Building 22. Several terrestrial targets are visible 
from this lab, including trees at various ranges, a water tower at ~5 km range, and 
power line towers near the local horizon, approximately 1 2 to 1 5 km range. Figure 42 
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shows an altimetry measurement to a water tower at ~5 km. System parameters for 
this measurement are listed in Table 1 2. This nighttime measurement was made with 
35 mW of average power and an integration time of one-second. A 4095 bit PN code was 
used at a bit rate of 1 MHz. The horizontal measurement resolution was 150 m and the 
data has not been smoothed nor range-corrected. The background count rate was 300 
counts/second during the measurement. 



0 5 10 15 

Range (km) 

Figure 42. Nighttime, horizontal path lidar measurements to a water tower with 35 
mW average transmit power and a 1 second integration time. 

Table 1 2 

Lidar system parameters for lidar measurement shown in Figure 42. 


Avg. laser power 

35 mW 

PN code length 

4095 bits 

Count rate 

220 Kcounts/sec 

Background count rate 

300 counts/sec 

Transmitter divergence 

< 1 00 urad 

Receiver FOV 

1 60 urad 

Bit rate 

1 MHz 

Integration time 

1 sec 
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The calculated signal-to-noise ratio (SNR) for the water tower return in Figure 42 is 


SNR = 



1.34 x IQ 5 counts 
V 2 .I x 10 5 counts 


= 292. 


( 99 ) 


Here S n is the photoelectron count of the correlation peak corresponding to the water 
tower and the denominator of Eq (99) is the sum of all counts in the histogram. Using 
the single scattering lidar equation Eq. (19) and the lidar system parameters, Table 12, 
the expected photoelectron count rate for a 5 km horizontal path at sea level from a 
target with a diffuse surface reflectivity of 0.5 is calculated to be N = 5 x 10 5 
counts / sec. A system transmission constant, = 0.5, and atmospheric attenuation 
coefficient K = 0. 1 km~ x are assumed. This attenuation coefficient is the sum of the 
aerosol scattering and absorption terms and a molecular scattering term for mid- 
latitude summer with 23 km visibility. The molecular absorption term has been 
assumed negligible. The total counts received in the one-second integration period was 
N = 2.1 x 10 5 counts. The observed count is approximately a factor of two lower than the 
calculated count rate. The lower rate is likely attributable to errors in the assumed 
values of the attenuation coefficient K, the system transmission and the diffuse 
reflectivity of the water tower. A better estimate of the target reflectivity is ~0.9 
which is based on the diffuse reflectivity of magnesium oxide paint on steel. 62 However 
this would lead to a larger discrepancy between observed and estimated count rates. 

Figure 43 shows a lidar measurement to a terrestrial target at approximately 1 3 
km. This measurement was made at night with 35 mW of average power and a one- 
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second integration time. A 255 bit PN code was used at a 1 MHz bit rate. The lidar 
system parameters for this measurement are listed in Table 1 3. 



Figure 43. Nighttime lidar measurements to a power-line tower with 35 mW average 
transmit power and a 1 second integration time. 


Table 1 3 

Lidar system parameters for lidar measurement shown in Figure 43. 


Avg. laser power 
PN code length 
Count rate 

Background count rate 
Transmitter divergence 
Receiver FOV 
Bit rate 

Integration time 


35 mW 
255 bits 

>150 Kcounts/sec 
5.2 Kcounts/sec 
< 1 00 urad 
320 urad 
1 MHz 
1 sec 
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The calculated SNR for the target return in Figure 43 is 


SNR = 



3.4 x 10 4 counts 


= 175. 


4 x 10 4 counts 


(100) 


The expected photoelectron count rate for this measurement as calculated by the lidar 
equation, using the system parameters in Table 13 and atmospheric parameters as in the 
previous calculation, is 1.54 x 10 4 counts / sec. The total received counts in the one- 
second measurement was 9.4 x 10 4 counts. The target is believed to be a power line 
tower with a large cross section central support. Daytime observations with binoculars 
reveal a series of these large power line towers in the target area, all painted with a 
highly reflective white paint. The only other targets in the area are trees. Previous 
lidar measurements to trees have never yielded such strong correlation peaks. The 
observed count rate is approximately a factor of six higher than the calculated count 
rate. It is possible that the strong lidar return is due to a high reflectivity paint used on 
power line towers to improve their visibility to low flying aircraft. The paint likely 
has characteristics which enhance its specular reflectivity. 

5.2 Cloud & aerosol lidar measurements 

The single-color, direct-current modulated, AIGaAs lidar system was also used to 
make lidar measurements of clouds and aerosols over slant paths. Once alignment was 
verified over a horizontal path, the system was pointed above the horizon, at several 
different elevation angles, and permitted to integrate for many minutes. Figure 44 
shows a measurement taken on the night of January 20, 1993 with 20 minutes of 
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integration time. The laser was operated at 35 mW average power and modulated with 
4095 bit code at a 1 MHz bit rate. The lidar system parameters are listed in Table 1 4. 
The data has been smoothed with a two-bin running average. Rayleigh scattering is 
evident up to approximately 9 km. Although that night was particularly clear with 
numerous stars visible, a thin layer of cirrus clouds is evident in the data at an altitude 
of approximately 1 3 km. This data has been range-corrected but not corrected for the 
overlap function. The calculated SNR at the peak of the cirrus return, prior to range- 
correction, is 


SNR = 


4095 

2 "< 

i = 0 


3.2 x 10 4 counts 
V5.1 x \0 6 counts 


= 13.9 . 


001 ) 
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Photoelectron Counts*10 6 *km 2 

Figure 44. Slant path lidar measurements to thin cirrus clouds at night. 

The Rayleigh backscatter, visible in Figure 44 to 9 km, may be compared with a model 

of the atmospheric molecular backscatter and the total attenuation coefficient at 

0.82 flm. Figure 45 shows a plot of the product of the volume backscatter coefficient and 

the attenuation coefficient. 63 
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Figure 45. Model of atmospheric backscatter and total attenuation coefficient vs 
altitude. [Spinhirne] 

The lidar data in Figure 44 does not show the increased boundary layer particulate 
scattering seen in the model. This is due to the lidar transmitter not completely 
overlapping the receiver until ~2 km altitude. However, the slopes of both the lidar data 
and the model atmosphere are similar in the region from 2-10 km. This comparison is 
shown in Figure 46. 


110 





Height (km) 

Figure 46. Comparison of model atmosphere with lidar data shown in Figure 44. 

Figure 47 shows a 20 minute lidar measurement taken immediately following the 
measurement shown in Figure 44. System parameters are listed in Table 14. A new 
cirrus layer forming below the layer in Figure 44 is apparent. 


Table 14 

Lidar system parameters for cirrus cloud lidar measurements. 


Avg. laser power 

35 mW 

PN code length 

4095 bits 

Count rate 

7-9 Kcounts/sec 

Background count rate 

1 40 counts/sec 

Transmitter divergence 

< 1 00 urad 

Receiver FOV 

1 60 urad 

Bit rate 

1 MHz 

Integration time 

1 200 sec 


Ill 


Lidar Signal 
(Photoelectron counts) 





Figure 47. Slant path lidar measurements to thin cirrus cloud layers at night. 

5.3 Integrated path water vapor DIAL measurements 

The lidar system depicted in Figure 37 has been used tQ make integrated path 
DIAL measurements of water vapor over a horizontal path through the boundary layer at 
night. Integrated path measurements return the average water vapor density over the 
measurement path, and allow greater measurement sensitivity at the expense of range 
information. For the on-line measurement, the AIGaAs laser was locked to an absorption 
line at 81 1 .61 7 nm. This line has a halfwidth, y , of 0.0837 cm -1 at STP and a line 

strength, 5, of 2.54x10 ~ 2A cm~' I molecules ■ cm~ 2 . The off-line measurement, made 
after the on-line, used the same laser actively stabilized at 81 1.640 nm, 23 pm away 
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from the absorption line. Figure 48 shows a 75 second integration for each on-line and 
off-line measurement. System parameters for the measurement are listed in Table 1 5. 
The largest peak in the figure is due to the reflection from a water tower at 5.1 km. 



Range (km) 

Figure 48. Integrated path DIAL measurements of water vapor over a horizontal path 
at night. Error bars are shown but very small due to strong signal. 


Table 1 5 

System parameters for water vapor lidar measurements shown in Figure 48. 


Avg. laser power 

25 mW 

PN code length 

1023 bits 

On-line count rate 

10.5 Kcounts/sec 

Off-line count rate 

23 Kcounts/sec 

Background count rate 

1 50 counts/sec 

Transmitter divergence 

<100 urad 

Receiver FOV 

1 60 urad 

Bit rate 

1 MHz 

Integration time 

75 sec 
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The ratio of the on-line to off-line peak is 0.268 and represents the integrated path 
attenuation due solely to water vapor. The integrated path absorption coefficient may be 
estimated using the Beer-Lambert attenuation law, 




= 0.268. 


( 102 ) 


This yields a measured absorption coefficient, k H20 , of 0.1 29 km- 1 . At the time of the 
measurement, the mean surface pressure reported by a local airport was 1025 mb with 
a vapor pressure of 5.64 mb. The absorption coefficient at line center for this vapor 
pressure, assuming a pressure broadened Lorentz line shape, is defined as 


k a -n -S 0 /K- Yo . 


003) 


and is calculated to be 0.1 38km- 1 . The measured absorption coefficient is 6.5% lower 
than calculated. The error is due to four possible sources ( 1 ) the lidar is separated 30 
km from the airport, where the reported surface and vapor pressure were measured, 

(2) the 0.8 mA sinusoidal current dither on the laser diode which causes a sinusoidal 
wavelength variation of ±1.3 pm around the absorption line center, (3) the frequency 
offset involved with manually closing the feedback loop, and (4) interference from a 
weak absorption line very close in wavelength to the absorption line at 81 1 .61 7 nm. 

The first error source may be minimized by making humidity measurements at 
the lidar’s location. The second error source, which leads to an underestimation of the 
water vapor density, may be accounted for by calculating the average or effective 
absorption coefficient rather than the line center absorption coefficient. The third 
error source, which also leads to an underestimation of the water vapor density, is more 
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difficult to account for in the data and is best dealt with by using the third derivative of 
the absorption line as the locking signal. Error sources (2) and (3) are dealt with in 
section 5.5, Error analysis. The fourth error source should lead to an overestimation of 
absorption and is best dealt with by selecting better isolated lines. 

5.4 Range resolved water vapor DIAL measurements 

The AIGaAs water vapor DIAL system was pointed above the horizon, to several 
different elevation angles, in order to eliminate returns from terrestrial targets. Again, 
both on-line and off-line data sets were accumulated with the laser wavelength 
stabilized throughout each one-minute measurement. Due to the proximity of the paths, 
the horizontal range-resolved humidity measurements should agree with the average 
humidity estimated from the integrated path measurements of the previous section. 
Figure 49 shows a single on-line measurement from the night of October 28, 1 993. The 
lidar was pointed 2-3° above the horizon. The AIGaAs laser was actively stabilized to the 
81 1.617 nm absorption line throughout the one-minute measurement. System 
parameters for the lidar measurements are listed in Table 16. 


Table 1 6 

System parameters for range-resolved DIAL measurements of atmospheric water vapor. 


Avg. laser power 20 mW 

PN code length 255 bits 

On-line count rate 0.38 Kcounts/sec 

Off-line count rate 0.55 Kcounts/sec 

Background count rate 50 counts/sec 

Transmitter divergence <100 urad 

Receiver FOV 1 60 urad 

Bit rate 1 MHz 

Integration time 60 sec 
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Figure 49. Range resolved lidar measurement of water vapor with laser tuned to line 
center of the 81 1 .61 7 nm water vapor absorption line. Trajectory is 2-3* above the 
horizon with a one minute integration time and 20 mW average power. 


This data has not been range corrected or smoothed. From the figure, it is evident that 
the transmitter beam is not completely overlapped by the receiver FOV until a range of 
nearly 1 .0 km. The exponential attenuation of the lidar signal from 1 to 3 km agrees 
with the 1 /R 2 losses and the Beer-Lambert attenuation law. The off-line measurement, 

Figure 50, was made shortly after the on-line measurement. The integration time was 
one-minute and the trajectory was 2-3° above the horizon. The laser was actively 
stabilized to 81 1.635 nm for the duration of the one minute measurement. 
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Figure 50. Range resolved lidar measurement with laser tuned to 81 1.635 nm, away 
from the water vapor. Trajectory is 2-3° above the horizon. 



This data has not been range corrected or smoothed. For ease of comparison, Figures 49 
and 50 are replotted on the same axes, Figure 51. The difference between the on-line 
and off-line lidar signals is primarily the absorption due to water vapor and to a lesser 
extent, the difference in transmitted power caused by tuning the laser wavelength with 
bias current. If we assume a slope efficiency of 0.3 mW/mA for the AIGaAs lidar 
transmitter, Figure 38, then the difference in on-line to off-line average transmitted 
power is ~1.5 mW or approximately a 5% change in the transmitted power. The 
difference in the peak photoelectron counts between on-line and off-line, at 1 km range, 
is -25%. Therefore, the change in transmitted power accounts for -20% of the 
difference in the on-line to off-line lidar signals. 
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Figure 51 . Overlay of one-minute on-line and off-line range resolved water vapor 
DIAL measurements. 


By increasing the integration times, the signal to noise ratio (SNR) of the 
measurements may be improved. To retain some temporal information however, 
individual one or two minute histograms are accumulated consecutively rather than as 
one long histogram. Afterwards, the histograms are added together and then cross- 
correlated with the PN code. Figure 52 shows on-line and off-line lidar profiles of two 
twelve minute measurements. The system parameters are listed in Table 1 7. The laser 
was actively stabilized to the 81 1.617 nm absorption line for the on-line measurement 
and to 81 1.640 nm, ~23 pm away from the absorption line center, for the off-line 
measurement. The correlation functions of twelve one minute histograms were summed 
together for each lidar profile. 
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Table 1 7 


System parameters for water vapor lidar measurements shown in Figure 52 & 53. 


Avg. laser power 

25 mW 

PN code length 

255 bits 

On-line count rate 

0.38 Kcounts/sec 

Off-line count rate 

0.55 Kcounts/sec 

Background count rate 

50 counts/sec 

Transmitter divergence 

<100 urad 

Receiver FOV 

1 60 urad 

Bit rate 

1 MHz 

Integration time 

60 sec 


Since the data accumulation (histogramming) and cross-correlation are linear 
operations, the correlation functions may be directly summed to yield the same result as 
adding the histograms prior to correlating. 
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Figure 52. On-line and off-line range resolved water vapor lidar measurements. 
Each lidar profile is the sum of 1 2 one minute integrations. 
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In Figure 53, the lidar data from Figure 52 has been range corrected and plotted 
on a logarithmic scale. The range correction involves multiplying the photoelectron 
count in each range bin by the square of the distance from the receiver. If the dominant 
signal loss factor of the off-line signal is the 1/ R 2 losses then the slope of the off-line 
signal should then be very nearly zero. Since the dominant attenuation of the on-line 
signal is due to 1 / R 2 losses and absorption by water vapor then the on-line signal should 
have a nearly constant negative slope, which is evident in the figure. 



Figure 53. Horizontal path water vapor DIAL data, from Figure 52, range corrected 
and plotted on a logarithmic scale. Integration time was 1 2 minutes. 


It should be noted that the error bars which previously were constant in size have also 
been multiplied by the square of the range. 
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Water vapor lidar data was taken on the night of November 4, 1 993 and is shown 


in Figure 54. This data has also been range corrected. System parameters for these 
measurements are listed in Table 1 8. 



Figure 54. Horizontal path water vapor DIAL data, range corrected and plotted on a 
logarithmic scale. Integration time was 1 1 minutes. 


Table 1 8 

System parameters for lidar measurements shown in Figure 54. 


Avg. laser power 

20 mW 

PN code length 

255 bits 

On-line count rate 

6.2 Kcounts/sec 

Off-line count rate 

9.0 Kcounts/sec 

Background count rate 

80 counts/sec 

Transmitter divergence 

<100 urad 

Receiver FOV 

1 60 urad 

Bit rate 

1 MHz 

Integration time 

1 1 minutes 
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The higher SNR in Figure 54, seen as smaller error bars, is due to an improvement in 
the receiver optical alignment between the data taken on October 28 and November 4, 

1 993. This improvement resulted in more than an order of magnitude increase in signal 
count rate. Figure 55 shows range-resolved water vapor DIAL measurements taken on 
the night of November 1 0, 1 993. Each on-line and off-line measurement was integrated 
for 20 minutes with the lidar pointed 2-3° above the horizon. The data has been range 
corrected. The system parameters for these measurements are listed in Table 1 9. 


Table 1 9 

System parameters for water vapor lidar measurements shown in Figures 55 & 56. 


Avg. laser power 25 mW 

PN code length 1 023 bits 

On-line count rate 7.4 Kcounts/sec 

Off-line count rate 10.0 Kcounts/sec 

Background count rate 1 50 counts/sec 

Transmitter divergence < 100 urad 

Receiver FOV 1 60 urad 

Bit rate 1 MHz 

Integration time 20 minutes 


A lidar measurement through the boundary layer should experience a nearly constant 
extinction coefficient. This should yield a constant slope when range corrected and 
plotted on a logarithmic scale, which is evident in the figure. 
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Figure 55. Horizontal path, range-resolved DIAL measurements of atmospheric water 
vapor made on the night of 1 1/10/93. Integration time was 20 minutes for each 
profile. 

The water vapor absorption coefficient was estimated from the slopes of the on-line and 
off-line measurements in Figure 54 to be 0.109 km~ x . At the time of the lidar 
measurement, the mean surface pressure was 1025 mb and the vapor pressure was 
5.41 mb. This vapor pressure corresponds to a water vapor number density of 
1. 37 x 10 17 molecules / cm 3 . The calculated absorption coefficient for the 811.617 nm 
absorption line, using this number density, is 0.132 km ~ x . The lidar measurement is 
1 7.4% lower than the calculated which is likely due to the same error sources listed in 
section 5.3. These error sources are addressed in section 5.5. 

Estimates of range-resolved water vapor number density can be made by 
applying the DIAL equation Eq. (6) to the data shown in Figure 55. However, the data is 
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somewhat noisy so a two-bin (300 m resolution) running average has been performed 
and the data plotted in Figure 56. 



Range (km) 

Figure 56. Water vapor DIAL data from 1 1/10/93 (Figure 55) smoothed in range 
with a two-bin running average. Range resolution = 300 m. 


The range resolved water vapor number density can be estimated by using the DIAL 
equation, 


n(AR) = 


1 

2AoA R 


J’ 


(104) 


where A <7 = 9.659 x 10 -24 cm 2 / molecules and A R = 300m. Applying Eq. ( 1 04) to the 
lidar data in Figure 56 yields an estimate of the range-resolved number density of water 
molecules in the atmospheric path. The zero and negative values of n(AR) prevent the 
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data from being plotted on a log axis. Since zero and negative values of atmospheric 
water vapor number density are not realistic they have been dropped and the remaining 
data plotted in Figure 57. Missing data has been interpolated. The mean water vapor 
density over the atmospheric path, calculated from all data points including negative and 
zero values, is 0.48 x W 1 molecules / cm 3 . This average water vapor density may be 
taken as the lower bound. The mean water vapor density calculated with only the 
positive values of n(AR) yields 1.49 xlO' 1 molecules / cm 3 , which may be taken as the 
upper bound 



Figure 57. Range-resolved water vapor number density estimated from lidar 
measurements made on 1 1/10/93. Negative and zero values of water vapor density have 
been dropped. Trace connecting data points interpolates between dropped data points. 
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Further smoothing of the lidar data should eliminate the zero and negative values of 
n(AR) at the expense of range information. As the number of bins added together 
increases, the estimated water vapor density should approach the mean absorption 
coefficient calculated by fitting a line to the same data. The mean absorption coefficient 
was estimated to be 0.109£m -1 . This absorption coefficient corresponds to a mean 
water vapor number density of 1.13 x 10 11 molecules I cm? or a mixing ratio of 
2.73g/kg. 

At the time of the water vapor DIAL measurement, the water vapor mixing ratio, 
measured at BWI Airport, was reported to be 3.47g / Agwhich corresponds to a number 
density of 1.37 x 10 molecules / cm 3 . The lidar underestimated the water vapor 
mixing ratio by 2 1 % when compared with the humidity measurements reported by BWI. 

5.5 Summary of water vapor DIAL measurements 

The results of the water vapor DIAL measurements have been tabulated and are 
listed in Table 20. Columns 2 & 3 present the water vapor absorption coefficient 
calculated from the mean surface pressure and vapor pressure as reported by Baltimore 
Washington International Airport and Andrews Air Force Base (ADW). Column 4 shows 
the water vapor absorption coefficient estimated from the AIGaAs lidar measurements. 
Column 5 shows the difference between the measured and calculated values of the water 
vapor absorption coefficient. The arrows indicate whether the lidar estimate was higher 
(ft) or lower (li) than the reported values. With only two exceptions, the AIGaAs lidar 
measured a lower water vapor absorption coefficient than was reported by BWI or ADW. 
The errors which would cause the AIGaAs lidar to underestimate the water vapor 
absorption coefficient are addressed in Section 5.6. 


126 



Table 20 


Summary of water vapor lidar measurements made from 1 0/28/93 to 11/1 0/93. The 
denotes intermittent rain showers during the measurements. 


Date 

BWI 

(km -1 ) 

ADW 

(km' 1 ) 

£=3SB "" l AlGaAs Lidar 
Estimate (km' 1 ) 

Difference(IT or IT) 
BWI/ ADW 

11/10/93 

0.145 

0.139 

0.129 integrated path 

11%/ 7% (IT) 


0.138 

0.133 

0.1 10 range resolved 

20% / 17% (IT) 

11/09/93 

0.143 

0.150 

0.156 range resolved 

9%/ 4% (IT) 


0.143 

0.170 

0.125 integrated path 

13%/ 26% (IT) 

11/08/93 

0.149 

0.132 

0.129 range resolved 

13%/ 2% (IT) 


0.138 

0.126 

0. 1 37 range resolved 

1% / 9% (lift) 


0.143 

0.132 

0. 1 1 5 range resolved 

19%/ 13% (IT) 

1 1/04/93* 

0.263 

0.245 

0. 1 67 range resolved 

36%/ 31% (IT) 

11/02/93 

0.150 

0.143 

0.079 range resolved 

47%/ 44% (IT) 

10/28/93 

0.158 

0.119 

0.113 range resolved 

28% / 5% (IT) 


5.6 Error analysis 

Both the integrated path and range resolved DIAL measurements made with the 
AIGaAs lidar underestimated the water vapor absorption coefficient and the water vapor 
number density. Since there is a significant distance between BWI and GSFC it is not 
unusual for there to be a discrepancy between the measured and estimated water vapor 
density. However, if the discrepancy were solely due to the spatial variation of water 
vapor in the atmosphere ( i.e . no instrument or measurement technique effects), one 
would expect the lidar measurements to return with equal likelihood both higher and 
lower absorption coefficients compared with the meteorological reports. Over the period 
of water vapor measurements, Oct. 28 through Nov. 10, the AIGaAs lidar consistently 
measured lower absorption coefficients than were reported at either Baltimore 
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Washington International Airport or Andrews Air Force Base. This implies a bias in the 
instrument or measurement technique which must be dealt with as an error source. 

There are two known error sources in the present DIAL technique which result in 
an underestimation of the absorption coefficient. They are ( 1 ) the manual zeroing of the 
lockin amplifier to the center of the absorption line and (2) the wavelength detuning 
effects of the dither current to the laser. 

The offset error induced by manually zeroing the lockin amplifier may be 
estimated by from Figures 61 and 62. This error may be represented as the ratio of 
short term lockin voltage fluctuations to the first derivative amplitude. From Figure 
62, the maximum change in lockin voltage between two consecutive data points is 
approximately 1 mV. The amplitude of the first derivative, from Figure 61, is ~ 17 mV. 
This ratio is 

= 0.058 (105) 

\lrnV 

which implies an offset error of 0.058 times the Doppler linewidth. The slope of line 
center of the first derivative curve, Figure 61 , is 


11 mV - 24. 3 mV / mA . 
0.7mA 


(106) 


The estimated offset error of ~1 mV pp corresponds to (lmV)/(24.3mV / / mA) = 0.04mA 
which translates to a wavelength fluctuation of 

(0.04mA)- (3.3pm/ mA) = 0.13 pm , (107) 
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or approximately one-tenth of a linewidth, which corresponds to less than 2% change in 
the absorption coefficient. 

The most significant error source is the wavelength detuning effect of the 0.8 mA 
dither current used to frequency lock the laser to the water vapor absorption line. By 
scanning the laser wavelength sinusoidally across the center of the absorption line, the 
laser experiences a lower effective absorption coefficient. This results in the AIGaAs 
lidar underestimating the water vapor absorption coefficient and, hence, the water vapor 
number density. 

By calculating this effective absorption coefficient, this error may be estimated 
and subsequently corrected for in the data. The effective, i.e. time averaged, absorption 
coefficient can be estimated by integrating the product of the Lorentz absorption line 
shape with the probability density of the sinusoidal wavelength dither signal. The 
Lorentz line shape is described by the equation, 



s r/2 
n { v ~ v o) 2 + {r/ 2 ) 2 ' 


(108) 


where S is absorption transition strength, y is the absorption linewidth, and v is the 
frequency. Figure 58 shows a plot of the normalized Lorentz linewidth with a halfwidth 
equal to unity and line center v 0 = 0. 
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Figure 58. Plot of Lorentz lineshape with halfwidth a = 1 . 


The probability density of a sine wave amplitude with peak amplitude 0.5 is given by® 4 


p(v) = 


1 Ivl < 0.5 

7T-^( 0.5) 2 -V 2 

0 |v|> 0.5 


(109) 


and is shown in Figure 59. 



dither amplitude - arb units 

Figure 59. Plot of probability density of sine function amplitude with amplitude = 0.5. 
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Numerical integration of the probability density function over the interval, 

-0.5 < v < 0.5, approaches unity as the limits approach ±0.5. The time averaged 
absorption coefficient cr(v) is the integral of the product of the normalized linewidth and 
probability density is represented as 


1 <+0.5 

1 f 

r 1 l. 

1 

n J 

>-0.5 

_v 2 + lj 

_Vo.25- v 2 _ 


\dv. 


010 ) 


For the above values, numerical integration of Eq (1 10) yields a value of 0.89. This 
implies that the absorption coefficient, estimated from the lidar data, will be 
underestimated by ~1 1% when the laser is dithered with an amplitude of one-half the 
linewidth. Figure 60 shows the normalized absorption coefficient as a function of dither 
amplitude. As expected, the absorption coefficient is unity for zero dither amplitude and 
decreases as the dither amplitude increases. As the dither amplitude grows to many 
times the linewidth, the absorption coefficient asymptotically approaches zero. 


131 



Wavelength Dither Amplitude 
(fraction of absorption linewidth) 


Figure 60. Normalized effective absorption coefficient vs wavelength dither 
amplitude. The peak dither wavelength is given as a fraction of absorption linewidth. 

The applied dither current was typically (40mV)/(50fl) = 0.8mA. This current 
corresponds to a peak wavelength dither of (0.8mA) -(3.3pm/ mA) = 2. 64 pm or 
±1.32pm . The Lorentz or pressure broadened linewidth of the 81 1 .61 7 nm absorption 
line is 10.98pm . This implies that the dither amplitude was approximately one- 
quarter of the linewidth. Calculating the integral in Eq (1 10) with a dither amplitude of 
0.25 yields a value for the normalized absorption coefficient of 0.97. This implies that 
the water vapor coefficients, and subsequently the water vapor number densities, 
estimated in sections 5.3 and 5.4 should be underestimated by ~3%. 
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Fluctuations in the wavelength of the laser represent another error source in the 
DIAL measurements. As the laser drifts away from the center of the absorption line, the 
laser will experience a reduced absorption coefficient which will result in an 
underestimation of the water vapor number density. The laser wavelength stability can 
be estimated by taking a ratio of the peak-to-peak excursions of the lockin error voltage 
and the peak-to-peak amplitude of the first derivative curve of the 811.617 nm 
absorption line. Figure 61 shows a first derivative scan of the absorption line at 
811.617 nm made with the absorption cell and AIGaAs lidar transmitter laser. The scan 
was made in an absorption cell with a 1 0 m path length which was evacuated and 
backfilled with ~26 Torr of water vapor. The first derivative curve of the water vapor 
absorption line is centered at 92.9 mA bias current and the linewidth is ~0.7 mA. This 
linewidth corresponds to 2.3 pm or 0.035 cm -1 (FWHM). The halfwidth, 0.017 cm -1 , 
agrees with the Doppler halfwidth calculated in Eq (92). 

The fluctuations surrounding the first derivative curve in Figure 61 may be 
attributed to an etalon effect created in the absorption cell. The curve in Figure 61 
represents the sum of the etalon effect and the first derivative absorption signal at 
81 1.617 nm. Although the laser is stabilized to the absorption line at 81 1.617 nm, 
Figure 62, the etalon effect could cause an offset error of the laser center wavelength 
which would lead to an error in estimating the water vapor absorption coefficient. 

The peak to peak excursion of the first derivative curve in Figure 61 is 1 7 mV. 

A scan of lockin error signal and laser wavelength, made during the water vapor DIAL 
measurements (Figure 55, data set 1 1 lOef), is shown in Figure 62. The top trace 
shows the laser wavelength during the lidar measurement. The bottom trace is the 
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lockin error voltage and shows a peak to peak excursion of 3 mV. The ratio of lockin 
error voltage to the first derivative amplitude is 


3mV 
17 mV 


0.17. 


(Ill) 



Figure 61. First derivative scan of 811 .61 7 nm water vapor absorption line. 

This ratio implies that the laser fluctuated less than one-fifth of the absorption 
linewidth during the DIAL measurement. The slope of line center of the first derivative 
curve, Figure 61 , is 24.3 mV / mA. The laser fluctuated ~3 mV pp during the water 
vapor DIAL measurement, which corresponds to ( 3mV)/(24.3mV / mA) = 0.123mA. 

This translates to a wavelength fluctuation of 

(0.123mA) (3.3/?m/mA) = 0.4 pm . (112) 
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Wavelength (nm) 




Therefore the laser center wavelength fluctuated 0Apm/2.3pm, or 0. 1 7 of the 
absorption linewidth, which is in agreement with the stability estimate from Eq (1 1 1). 
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Figure 62. Frequency stability of AIGaAs laser during on-line water vapor lidar 
measurement of 11/1 0/93. Top trace is wavelength and bottom trace is lockin error 
voltage. 


It is important to note that the Doppler width of the absorption line at 81 1 .61 7 
nm is approximately one-fifth of the Lorentz or pressure broadened linewidth, 

(O.Qlll cm' 1 }/ (0.0837 cm~ 1 } = 0.21 . The laser center wavelength was stable to better 
than one-fifth of the Doppler width, therefore, with respect to the Lorentz linewidth, 
the laser was stable to (0.2) (0.2) or 0.04 times the Lorentz linewidth. The absorption 
coefficient, Eq (108), can be calculated for a laser frequency drift from line center of 
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Wavelength (nm) 




0.04, i.e. (v - v c ) 2 = 0.0016. The resultant absorption coefficient is 0.998 the value at 
line center, implying a DIAL measurement error of approximately 0. 1 6%. 


A summary of error sources and the estimated error in the absorption coefficient 
is presented in Table 2 1 . 


Table 21 

Summary of error sources in AIGaAs water vapor DIAL measurements and magnitudes of 
these errors in estimates of the water vapor absorption coefficient. 


Error Source 

Estimated error to 
absorption coefficient 

Physical Separation of lidar and ground 
truth humidity measurements 

-0-40% 

Dither current induced wavelength 
variation ~0.3 linewidth 

-3-5% 

Manual zeroing of lockin amplifier 
-0.25 linewidth from linecenter 

<2% 

Drift of laser center wavelength 
<0.04 of Lorentz linewidth 

«1% 

Interference from neighboring lines 

-0.2% 
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6. SUMMARY AND CONCLUSIONS 


Single color lidar measurements over horizontal paths to terrestrial targets 
were made at night. Horizontal path measurements to a water tower at 5 km were used 
to check and optimize the alignment of the lidar transmitter with the receiver and to 
compare the operation of the lidar system with theory via the single scattering lidar 
equation. Calculated performance agreed with experimental to within a factor of two 
which is reasonable considering the accuracy of assumed values for the target 
reflectivity and the total atmospheric attenuation coefficient. 

Nighttime cloud and aerosol lidar measurements were made over near vertical 
paths. Rayleigh backscatter was visible in 20 minute lidar measurements to an altitude 
of approximately 9 km. Thin cirrus clouds at 10 and 13 km altitudes were profiled with 
an average laser power of 35 mW and an integration time of 20 minutes. 

The cloud and aerosol lidar measurements made in these experiments represent a 
significant increase in measurement range and accuracy compared to previous 
measurements reported with PN code AIGaAs lidar systems. 

DIAL measurements of atmospheric water vapor have been made for the first 
time using a PN code modulated AIGaAs laser. An external EO modulator was used to 
impress the PN code onto the laser diode beam. Integrated path measurements were made 
to a terrestrial target over a one-way path length of ~5 km. The estimated water vapor 
absorption coefficient agreed to within 6.5% of the measured humidity reported by 
Baltimore Washington International Airport at the time of the lidar measurement. 
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Range resolved DIAL measurements of water vapor were made over a 4 km horizontal 
path. The estimated mean water vapor coefficient from these measurements agreed to 
within 20% of the measured humidity. The humidity values used to compare with both 
integrated path and range resolved water vapor lidar estimates were removed from the 
lidar system by approximately 30 km. 

Error sources in the DIAL measurements were identified and quantified where 
possible. The 30 km separation of the lidar system and meteorological station 
represents an indeterminate error. This error source may be reduced by making 
humidity measurements at the location of the lidar. The wavelength stability of the 
AIGaAs laser was measured and the resulting DIAL measurement error due to laser drift 
was estimated to be less the 0.2%. The most significant error source was the dither 
current used to frequency lock the AIGaAs laser to the water vapor absorption line. This 
error was a function of the dither amplitude. The resulting error ranged from 3% for a 
dither amplitude equal to one-quarter of the absorption linewidth to 60% for a dither 
amplitude equal to the full linewidth. However, this was a systematic error which, once 
quantified can be removed from the data by applying a multiplicative correction factor. 

The DIAL measurements of water vapor are the first reported measurements of 
atmospheric water vapor using a compact lidar based on PN code modulated AIGaAs 
lasers. 


6.1 Future work 

Work is underway to replace the AIGaAs laser transmitter with a 1 Watt CW 
master oscillator-power amplifier (MOPA) device. This could potentially eliminate the 



need for the external EO modulator by permitting current modulation of the power 
amplifier section of the MOPA device. This depends on whether modulating the power 
amplifier effects the spectral properties, i.e. linewidth and center frequency, of the 
master oscillator. The master oscillator would be frequency locked to a water vapor 
absorption line. 

The optical bandpass filter currently used in the lidar receiver (810 nm center 
wavelength, 5 nm FWHM) will be replaced with a custom optical bandpass filter with 
<0.1 nm bandwidth. This should permit daytime operation of the AIGaAs lidar. 

An improved histogramming circuit has been designed and will be implemented on 
a custom printed circuit board. This new histogramming circuit will eliminate 
measurement errors due to the existing histogrammer. 

Improved opto-mechanical design of the lidar system should result in a more 
compact lidar and better alignment stability. Also, improvements to the thermal design 
of the laser diode header are needed to reduce temperature induced fluctuations of the 
laser diode wavelength. This should also improve the frequency locking of the AIGaAs 
laser to water vapor lines. 

An improved frequency locking algorithm is being developed. This will 
incorporate a true proportional-integral-derivative (RID) feedback loop which should 
improve both the short term, < 1 minute, and long term, ~1 hour, frequency stability. 
Also, techniques will be investigated to permit simultaneous on-line and off-line DIAL 
measurements. It is not clear at this time whether a single laser, frequency hopping 
between the on-line and off-line wavelengths, is simpler than two separate laser 
transmitters, orthogonally polarized, operating at the on-line and off-line wavelengths 
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simultaneously. This latter technique would require two photon count detectors and a 
polarization beam splitter in the receiver path. 

The multipass absorption cell will be investigated and the etalon effect minimized 
or eliminated. The frequency locking technique used to stabilize the laser to water vapor 
absorption lines is susceptible to this etalon effect. Improved frequency locking to more 
lines, and weaker lines, should be realized with elimination of the etalon effect. 

A prototype AIGaAs altimeter will be assembled and tested. The AIGaAs altimeter 
will use a single element, 150 mW CW AIGaAs laser modulated with a 2047 bit PN code 
at >100 MHz. The range resolution will be less than 1.5 meters. This prototype 
altimeter could be flown on a NASA aircraft and compared with existing laser altimeter 
instruments. 
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